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Appendices 


INTRODUCTION 

Adjuvant  chemotherapy  reduces  the  risk  of  breast  cancer  death  for  essentially  all 
stages  of  the  disease,  although  the  absolute  benefits  of  chemotherapy  are  small.  The 
development  of  many  new  effective  agents  may  provide  several  improved  options  for 
many  women.  A  major  challenge  in  breast  cancer  treatment  is  to  determine  which 
chemotherapeutic  agent  will  provide  the  most  benefit  for  an  individual  patient.  Effective 
chemotherapy  inhibits  cell  proliferation  and  induces  apoptosis.  In  previous  studies  of 
cells,  it  has  shown  that  elevated  levels  of  choline-containing  compounds  (tCho)  are 
associated  with  proliferation  (7),  and  that  changes  in  cellular  tCho  levels  occur  during 
malignant  transformation  and  progression  (2,  3)  and  in  altered  resistance  to  chemotherapy 
(4).  In  vivo  tumor  concentrations  of  tCho,  [tCho],  can  be  measured  non-invasively  using 
in  vivo  magnetic  resonance  spectroscopy  (MRS)  and  have  been  shown  to  be  linked  to 
malignant  proliferation  in  brain  tumors  (5).  The  objective  of  this  research  is  to  determine 
whether  in  vivo  MRS  measurements  of  [tCho]  can  identify  chemotherapy-induced 
growth  arrest,  and  possibly  apoptosis,  in  patients  with  locally-advanced  breast  cancer 
(LABC)  who  receive  primary  systemic  therapy  (PST),  also  known  as  neoadjuvant 
chemotherapy.  The  goal  is  to  determine  whether  changes  in  [tCho]  can  identify 
responders  and  non-responders  very  early  in  the  course  of  treatment  (within  24  hr  -  3 
wks).  The  study  employs  a  MR  scanner  having  a  high  magnetic  field  (4  Tesla)  which 
improves  the  sensitivity  for  [tCho]  measurement.  To  measure  tumor  size  changes  and  to 
locate  the  tumors  for  [tCho]  measurements,  the  protocol  also  includes  contrast-enhanced 
MRI  scanning.  MRS  measurements  of  [tCho]  are  compared  with  clinical  response  based 
on  tumor  size  changes  by  MRI  and  by  pathological  evaluation  of  the  surgically-resected 
tissue  after  PST.  This  research  is  expected  to  show  that  MRS  can  be  used  to 
individualize  and  optimize  adjuvant  chemotherapy  regimens  for  the  treatment  of  breast 
cancer. 


BODY 

In  the  period  of  time  since  gaining  approval  by  the  HSRRB  to  use  human  subjects 
(January  2002),  we  enrolled  35  patients  in  the  study  and  a  total  179  MRS/MRI  exams 
were  performed  on  these  subjects  for  the  purpose  of  accomplishing  Task  1:  Monitor 
changes  in  [tCho]  by  MRS  and  tumor  size  by  MRI  in  response  to  PST. 

After  completing  data  collection  of  the  first  13  subjects  who  completed  the  full 
course  of  PST,  we  analyzed  the  data  according  to  Task  2:  Determine  how  changes  in 
MRS  correlate  with  tumor  response  detected  by  MRI  and  with  pathology  findings.  The 
results  are  extremely  exciting,  since  the  MRS  measured  changes  in  [tCho]  appear  to 
predict  clinical  response  consistently  very  early  in  the  course  of  treatment.  Because  of 
the  substantial  correlation  between  changes  in  [tCho]  with  clinical  response,  we  decided 
to  report  these  novel  findings  before  completing  the  planned  accrual  (see  below).  These 
findings  are  to  be  published  in  November  2004  (b)  in  an  article  entitled  “Predicting 
Response  to  Neoadjuvant  Chemotherapy  of  Locally  Advanced  Breast  Cancer  with  In 
Vivo  'H  MRS:  A  Pilot  Study  at  4  Tesla”.  From  the  preprint  of  this  article  that  is  included 
in  the  Appendix,  it  can  be  seen  that  plots  in  figures  1,  3,  and  5  show  the  most  relevant 
data.  As  described  in  the  article  (b),  changes  in  [tCho]  between  the  baseline  scan 
(acquired  within  1  wk  before  starting  PST)  and  24  hours  after  the  first  cycle  of  PST 
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showed  a  significant  positive  correlation  with  change  in  tumor  size  (longest  dimension, 
LD)  after  4  cycles  of  doxorubicin-based  PST  (R=0.79,  p<0.001).  Opposite  trends  in 
[tCho]  changes  between  objective  responders  and  nonresponders  (p=0.0101)  were 
observed  already  by  24  hours  after  starting  PST. 

The  total  number  of  subjects  studied  (35)  was  less  than  the  target  number  (77). 
We  initially  powered  the  study  to  detect  changes  in  [tCho]  sufficient  to  distinguish 
between  patients  achieving  a  pathological  complete  response  from  those  with  residual 
tumor  left  at  definitive  surgery.  However,  with  the  high  accuracy  of  [tCho]  to  identify 
early  response  to  PST,  the  projected  sample  size  was  not  needed  to  accomplish  the  goals 
of  this  research.  As  described  above,  data  obtained  on  the  first  13  subjects  already 
provided  strong  evidence  of  ability  of  MRS  to  identify  early  response  to  PST  in  patients 
with  LABC.  The  presentation  given  at  last  year’s  RSNA  meeting  generated  much 
interest  among  radiologists  and  oncologists  in  the  field.  We  are  currently  getting 
inquiries  about  our  study  from  breast  cancer  imagers  and  researchers  from  multiple 
institutions  across  the  country.  In  late  Summer  2004,  a  teleconference  was  organized  to 
discuss  the  possibility  of  conducting  a  large  multi-center  clinical  trial  based  on  our  pilot 
study.  As  a  result,  an  initial  application  has  recently  been  submitted  to  ACRIN 
(American  College  of  Radiology  Imaging  Network),  an  agency  that  funds  these  types  of 
large  multi-center  trials.  We  hope  to  be  successful  in  obtaining  funds  so  that  the  multi¬ 
center  trial  can  begin  before  the  end  of  2005. 

The  MRS  methodology  used  to  quantify  tCho  in  breast  lesions  is  crucial  to  this 
project.  Although  the  development  of  this  methodology  was  not  one  of  the  specific 
Tasks  in  this  project,  the  data  generated  as  part  of  Task  1  could  be  used  to  help  validate 
our  tCho  quantification  procedure.  An  article  describing  our  quantitative  MRS 
methodology  entitled  “In  Vivo  Quantification  of  Choline  Compounds  in  the  Breast  with 
1H  MR  Spectroscopy”  was  published  in  2003  (7)  (see  Appendix).  In  addition,  in  the 
process  of  conducting  this  research,  we  solved  technical  problems  related  to  frequency 
shifts  that  previously  degraded  the  quality  of  MR  spectra  of  breast.  This  article  is 
currently  in  press  ( 8 )  (see  Appendix). 


KEY  RESEARCH  ACCOMPLISHMENTS 

•  We  have  shown  that  the  concentration  change  of  choline-containing  compounds 
in  the  tumors  (baseline  versus  24  hours  after  the  1st  cycle  of  chemotherapy) 
predicts  clinical  response  (as  determined  by  the  change  in  tumor  size  measured  by 
MRI;  baseline  versus  4  cycles  of  chemotherapy). 

•  These  results  suggest  that  MRS  can  be  used  clinically  to  distinguish  responders 
versus  nonresponders  early  in  the  course  of  neoadjuvant  chemotherapy.. 


REPORTABLE  OUTCOMES 

Manuscripts: 

S.  Meisamy,  P.J.  Bolan,  E.H.  Baker,  E.  Gulbahce,  L.I.  Everson,  M.T.  Nelson,  T.  Emory, 

T.  Tuttle,  D.  Yee,  and  M.  Garwood,  “Predicting  Response  to  Neoadjuvant  Chemotherapy 
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of  the  Locally  Advanced  Breast  Cancer  with  In  Vivo  1H  MRS:  A  Pilot  Study  at  4 
Tesla”,  Radiology,  in  press,  2004. 

P.J.  Bolan,  S.  Meisamy,  E.H.  Baker,  J.  Lin,  T.  Emory,  M.  Nelson,  L.I.  Everson,  D.  Yee, 
and  M.  Garwood,  “In  Vivo  Quantification  of  Choline  Compounds  in  the  Breast  with  1H 
MR  Spectroscopy”,  Magn.  Reson.  Med.,  50, 1 134-1 143, 2003. 

Measurement  and  Correction  of  Respiration-Induced  Bo  Variations  in  Breast  lH  MRS  at 
4  Tesla,  P.J.  Bolan,  P.-G.  Henry,  E.H.  Baker,  S.  Meisamy,  and  M.  Garwood,  Magn. 
Reson.  Med.,  in  press,  2004. 
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CONCLUSIONS 

Clinical  response  can  be  predicted  as  early  as  24  hours  after  the  1st  cycle  of 
doxorubicin-based  chemotherapy  using  in  vivo  MRS  measurement  of  choline-containing 
compounds  in  the  tumor.  These  findings  are  highly  significant  because  they  show  that 
MRS  can  distinguish  responders  versus  nonresponders  early  in  the  course  of  treatment 
and  suggest  the  use  of  MRS  clinically  to  individualize  treatments  for  maximizing  benefit 
and  to  rapidly  assess  efficacy  of  new  drugs. 
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Measurement  and  Correction  of  Respiration-Induced  B0 
Variations  in  Breast  aH  MRS  at  4  Tesla 

Patrick  J.  Bolan,*  Pierre-gilles  Henry,  Eva  H.  Baker,  Sina  Meisamy,  and 
Michael  Garwood 


Respiratory  motion  is  well  known  to  cause  artifacts  in  magnetic 
resonance  spectroscopy  (MRS).  In  MRS  of  the  breast,  the  dom¬ 
inant  artifact  is  not  due  to  motion  of  the  breast  itself,  but  rather 
it  is  produced  by  B0  field  distortions  associated  with  respiratory 
motion  of  tissues  In  the  chest  and  abdomen.  This  susceptibility 
artifact  has  been  reported  to  occur  in  the  brain,  but  it  is  more 
apparent  in  the  breast  due  to  the  anatomic  proximity  of  the 
lungs.  In  the  breast,  these  B0  distortions  cause  shot-to-shot 
frequency  shifts,  which  vary  an  average  of  24  Hz  during  a  typical 
*H  MRS  scan  at  4  T.  This  variation  can  be  corrected  retrospec¬ 
tively  by  frequency  shifting  individual  spectra  prior  to  averaging. 
If  not  corrected,  these  shifts  reduce  spectral  resolution  and 
increase  peak  fitting  errors.  This  work  demonstrates  the  arti¬ 
fact,  describes  a  method  for  correcting  it,  and  evaluates  its 
impact  on  quantitative  spectroscopy.  When  the  artifact  is  not 
corrected,  quantification  errors  increase  by  an  average  of  28%, 
which  dramatically  impacts  the  ability  to  measure  metabolite 
resonances  at  low  slgnal-to-noise  ratios.  Magn  Reson  Med 
52:000-000,  2004.  ©  2004  Wiley-Liss,  Inc. 

Key  words:  respiration;  susceptibility;  MRS;  breast  cancer  ? 

INTRODUCTION  r  Y  If 

Physiologic  motion  associated  with  respiratory  and  car¬ 
diac  cycles  is  known  to  cause  artifacts  in  magnetic  reso¬ 
nance  spectroscopy  (MRS).  The  situation  isl|nost  severe 
when  the  tissue  being  studied  experiences  g^oss  motibri, 
causing  the  voxel(s)  to  be  displaced  from  the  intended 
position.  This  is  common  in  thoracic  and  abdominal  MRSj 
where  techniques  such  as  triggering  (1,2),  breath-holding 
(3),  and  retrospective  data  rejection  (4)  can  be  used  to 
reduce  the  impact  of  gross  d||placement.  Whep  motions 
are  small  relative  to  the  voxel  sfze,  this  4i®laheip|rit  arti¬ 
fact  is  negligible.  However,  Ihefe  small-scale  tfotimilcail ' 
produce  phase  variations  \mwem  Jlqilfeei 

free  induction  decays  (FIDs).  If  these  are  not  corrected 
prior  to  averaging  the  FIDs,  destructive  interference  will 
cause  a  decrease  of  signal  intensities.  These  phase  varia¬ 
tions  can  be  corrected  by  retrospectively  phasing  individ¬ 
ual  spectra  prior  to  averaging.  The  reference  for  the  phase 
correction  can  be  derived  from  physiologic  monitoring  (5), 
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from  a  navigator  echo  (6,7),  or  from  resonances  in  the  data 
itself  (8-12).  Alternatively,  phase  variations  can  be 
avoided  altogether  by  using  physiologic  triggering  (13). 

MRS  has  been  increasingly  used  in  adjunct  with  breast 
MRI  for  diagnosing  cancer  and  monitoring  response  to 
cancer  treatments.  Malignant  breast  lesions  have  been 
shown  to  contain  elevated  levels  of  choline-containing 
compounds,  which  appear  as  a  single  resonance  at 
3.2  ppm  in  *H  MRS  (14).  Typical  breast  MRS  studies  are 
performed  with  the  subject  lying  prone  in  the  magnet  and 
with  a  coil  designed  to  support  the  breast.  Consequently, 
motion  of  the  breast  tissue  is  restricted,  and  gross  displace¬ 
ments  are  usually  less  than  1-2  mm.  This  motion  produces 
the  small  edge  artifacts  that  are  commonly  seen  in  subtrac¬ 
tion  images  of  dynamic  contrast  enhanced  breast  studies. 
Although  these  displacements  are  small  compared  to  the 
typical  *H  MRS  voxel  size  ( — 10  mm),  they  can  still  pro¬ 
duce  substantial  phase  variations. 

jAnother  effect  of  respiratory  motion  can  occur  even  when 
the  tissue  within  a  voxel  is  completely  immobile.  As  the 
li&gs  and  diaphragm  move  during  respiration,  they  create 
thhfe-dependint  Bp  variations  that  extend  far  from  the  chest 
(15,16).  This  was  demonstrated  by  Raj  et  al.,  who  showed 
that  respiratory  ihbtion  can  produce  significant  signal  fluc¬ 
tuations  in  EPI  iihages  acquired  in  a  phantom  placed  adja¬ 
cent  to  £  subject's  head  (17).  These  B0  fluctuations  produce 
frequency  shifts  in  acquired  spectra.  Previous  studies  in 
tein  MRS  have  shown  that  this  artifact  degrades  the  quality 
of  spectra  averaged  over  multiple  acquisitions  by  reducing 
-Spectral  resolution  and  increasing  quantitative  fitting  errors 
(ll,12,l|)Hn,the brain,  the  magnitude^ these  shifts  is  on  the 
order  opo^ppM^ffiippQ^MO:8-2  Hz  at  3  T  (12,18)  and 
fj.5— ■ 4  Hz  at  7  ®10||  Thl^shlfes jjj|re  small  relative  to  the 
Spectral  linewiflth,  s8r4ne  ^StraPdegradation  caused  by 
frequency  shifts  is  a  minor  effect  compared  to  the  signal  loss 
caused  by  phase  variations  (11).  However,  the  magnitude  of 
respiratory-induced  B0  shifts  increases  with  proximity  to  the 
chest  (16,19).  As  is  demonstrated  in  this  work,  the  magnitude 
of  these  shifts  in  the  breast  is  typically  0.1  ppm — 10  X  greater 
than  that  measured  in  the  brain. 

This  article  aims  to  measure  this  artifact  and  evaluate  its 
impact  on  quantitative  breast  MRS.  Breath-hold  imaging 
and  physiologic  monitoring  were  used  to  determine  the 
relationship  between  MRS  frequency  shifts  and  respira¬ 
tion.  Spectral  fitting  and  quantification  were  performed 
both  with  and  without  retrospective  frequency  correction 
to  demonstrate  that  fitting  errors  decrease  significantly 
when  frequency  correction  is  performed. 

METHODS 

All  measurements  were  performed  with  a  research  4-T 
system,  consisting  of  a  90-cm  bore  magnet  (model  4  T-900, 
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Oxford  Magnet  Technology,  Oxfordshire,  UK)  with  a  clin¬ 
ical  gradient  system  (Model  Sonata,  Siemens,  Erlangen, 
Germany)  interfaced  with  an  imaging  spectrometer  (Model 
Unity  Inova,  Varian,  Palo  Alto,  CA).  Several  different  sin¬ 
gle-breast  quadrature  transmit/receive  RF  surface  coils  of 
similar  design  were  used  to  accommodate  different  breast 
sizes  (20).  The  coils  were  mounted  onto  a  custom-built 
patient  table  designed  for  unilateral,  prone  breast  studies. 
In  several  examinations,  a  pneumatic  belt  was  wrapped 
around  the  thorax  just  below  the  breasts  to  acquire  chest 
expansion  data.  The  subject  population  for  this  study  in¬ 
cluded  both  normal  volunteers  and  patients  with  either  a 
suspicious  breast  lesion  or  a  known  breast  cancer.  The 
study  was  approved  by  our  institutional  review  board,  and 
informed  consent  was  obtained  prior  to  performing  each 
experiment.  All  processing  of  images  and  spectra  was  per¬ 
formed  using  Matlab  (The  Mathworks,  Natick,  MA). 

Assessment  of  Respiratory  Effects 

Breath-hold  imaging  was  performed  in  three  normal  sub¬ 
jects  to  evaluate  the  spatial  variation  of  respiration  arti¬ 
facts.  Single-slice  gradient  echo  images  (matrix  =  256  X 
128,  slice  thickness  =  5  mm,  field  of  view  =  14  cm, 
TR/TE  =  50/15  ms)  were  acquired  during  breath  holds  at 
both  maximum  inspiration  and  maximum  expiration, 
Maps  of  the  phase  difference  (Ao)  between  inspiration  and 
expiration  were  created  by  dividing  the  two  comfilex  im-  * 
ages  and  extracting  the  phase.  A  mask  for  the  high  signal- 
to-noise  region  was  created  by  thresholding  the  magnitude 
image  and  applying  spatial  smoothing  (using  Matllb’s  bv&i 
morph  function).  The  Ao  map  was  smoothed?  (mechan  fil¬ 
ter,  8X8  kernel)  and  manually  unwrapped  to  correct 
phase  discontinuities  greater  than  tt.  The  Ao  ihap  was  then  $ 
converted  into  a  frequency-shift  map  Av  ='|A6/(27rTE), 
expressed  in  hertz. 

Spectroscopy  Acquisition 

In  vivo  spectra  were  recorded  as  part  of  our  standard  breast 
MR  protocol,  consisting  of  a-dynamic  coijfrajj-enhanced 
MRI  study  followed  by  siitfjjle|yox§l  sl&Sfoscd^  iffi& 
localized  shimming  and  pjqw'CT^c 

spectra  were  acquired  usifig  loclhzanon  by  afiatiaticse- 
lective  refocusing  (LASER)  (21).  Water  suppression  was 
performed  using  variable  pulse  power  and  optimized  re¬ 
laxation  delays  (VAPOR)  (22).  TE  averaging  (23)  was  used 
to  reduce  lipid  sideband  artifacts  (TE  =  45-196  ms  in  N  = 
64  or  128  increments,  TR  =  3  s).  Each  FID  was  individu¬ 
ally  saved — no  averaging  was  performed  until  processing. 

A  fully  relaxed,  single-shot  spectrum  with  no  water  sup¬ 
pression  was  also  acquired  from  each  voxel  to  provide  a 
reference  water  resonance. 

Postprocessing 

Spectra  were  postprocessed  and  the  levels  of  total  choline- 
containing  compounds  (tCho)  were  calculated  using  a 
fully  automated  method  described  previously  (14)  and 
summarized  briefly  here.  For  postprocessing,  each  indi¬ 
vidually  saved  FID  was  automatically  corrected  for  DC 
offset  and  variation  of  zero-order  phase.  Frequency  shifts 
between  individual  FIDs  were  calculated  using  a  frequency- 


domain  cross-correlation  method  (18).  Each  FID  was 
shifted  in  frequency  by  multiplying  the  time-domain  data 
by  a  linear  phase  function,  exp(/-2iT-Av-f).  The  frequency 
shift  magnitude  Av  was  chosen  to  maximize  the  cross¬ 
correlation  (evaluated  in  modulus  mode  and  in  the  fre¬ 
quency  domain)  between  the  spectrum  and  a  reference 
spectrum  (arbitrarily,  the  first  in  the  series).  The  variability 
of  Av  over  all  FEDs  in  a  single-voxel  acquisition  was  char¬ 
acterized  by  its  range  R(Av)  and  SD  cr(Av). 

The  quantification  method  used  the  unsuppressed  water 
resonance  as  an  internal  reference  to  normalize  the  ampli¬ 
tude  of  the  tCho  resonance  and  obtain  an  estimate  of  the 
concentration  (denoted  [tCho])  expressed  in  molal  units 
(mmol /kg- water).  Spectral  fitting  was  performed  one  reso¬ 
nance  at  a  time  by  fitting  a  Voigt  line  shape  to  a  narrow 
band  of  the  frequency  domain  (0.4  ppm)  centered  around 
the  resonance.  The  fitting  error  (i.e.,  the  measurement  un¬ 
certainty)  was  estimated  using  the  Cram6r-Rao  minimum 
variance  bound,  which  was  calculated  using  the  covari¬ 
ance  matrix  provided  by  the  fitting  routine.  The  SD  of  the 
fitting  error,  <rtCho,  was  normalized  to  the  amplitude  of  the 
tCho  resonance  and  expressed  as  a  percentage.  This  error 
was  also  used  to  establish  the  detectability  threshold  as 
follows:  if  trtcho  was  greater  than  100%  then  the  tCho 
resonance  was  considered  undetectable. 

Quantitati  veSpectroscopy 

A:  retrospective  analysis  was  performed  to  evaluate 
Whether  the  frequency  correction  impacts  the  results  of 
spectral  quairtififcation.  In  vivo  breast  MRS  data  acquired 
fern  subjects  over  a  2-year  period  were  included  in  this 
;  sfijdy.  The  subjects  included  21  patients  with  biopsy-con- 
firniM  cdhcer  undergoing  neoadjuvant  chemotherapy,  46 
patiehts  with  a  suspicious  breast  lesion  evaluated  imme¬ 
diately  pre-  or  Jpostbiopsy,  and  5  normal  volunteers.  A 
total  of  715  TE-averaged  spectra  were  acquired  from  all  72 
subjects  in  140  studies.  A  subset  of  280  spectra  were  cho- 
,sen  for  this  retrospective  analysis  by  removing  9  spectra 
with  visible  artifacts,  151  spectra  with  lipid  content 
yggater  preyiqusly^^gcj^d  threshold  (14),  and 

j§ WB  spe|^'Withfac^  resonance.  These  280 

|lpectraliere  pi||ces^^phd%]i|Pt^ed  both  with  and  with¬ 
out  frequency  correction.  Omission  of  the  frequency  cor¬ 
rection  step  is  equivalent  to  the  common  practice  of  auto¬ 
matically  averaging  FIDs  prior  to  postprocessing. 

RESULTS 

Assessment  of  Respiratory  Effects 

Figure  1  is  a  representative  data  set  showing  sagittal  im¬ 
ages  from  one  breath-hold  study  of  a  normal  breast.  Mag¬ 
nitude  images  acquired  during  maximum  inspiration  (Fig, 
la)  and  maximum  expiration  (Fig.  lb)  are  very  similar.  The 
magnitude  difference  image  (Fig.  lc)  shows  that  the  gross 
motion  of  the  breast  was  small  during  respiration,  causing 
only  small  edge  artifacts.  The  phase  difference  between 
inspiration  and  expiration  was  converted  into  a  map  of 
frequency  shift,  Av,  shown  in  Fig.  1  day.  The  Av  map 
shows  a  strong  spatial  dependence,  with  the  greatest  effect 
(—70  Hz)  closest  to  the  chest  wall.  A  voxel  placed  in  the 
center  of  the  breast  would  have  experienced  a  shift  of 
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FIG.  1.  Demonstration  of  respiration  artifacts  by 
MR  imaging.  The  top  row  shows  sagittal  gradient- 
echo  images  of  a  subject’s  left  breast  during 
breath-holds  at  (a)  maximum  inspiration  and  (b) 
maximum  expiration.  The  magnitude  difference 
between  these  two  images  is  shown  in  c,  scaled 
by  a  factor  of  10,  showing  that  there  was  little 
gross  motion  of  the  breast  during  respiration.  The 
phase  difference  between  a  and  b  was  converted 
Into  a  Av  map,  shown  in  d,  with  contour  lines  every 
10  Hz.  The  frequency  shift  varied  from  <10  Hz  at 
the  skin  to  70  Hz  nearest  the  chest  wall. 


Av  (Hz) 


I  % 

10-20  Hz  between  maximum  inspiration  anil  expiration. 
Similar  Av  maps  were  reproduced  in  each^breath-jiold  , 
subject.  The  largest  frequency  shifts  were  neaiftst  the  cfeest 
wall  in  both  left  and  right  breasts  and  in  axial'|nd  sagittal 


base  provides  a  clear  example  of  the  impact  of  the  fre¬ 
quency  correction  method. 

Several  representative  examples  of  TE-averaged  acqui¬ 
sitions  "Sre  shown  in  Figs.  4b- d,  with  and  without  fre- 


orientations.  B 

The  relationship  between  MRS  frequency  shifts  and  thel 
respiratory  cycle  is  shown  in  Fig.  2.  In  the  figure,  the  dots 
indicate  frequency  shifts  measured  from  NMR  spectra  us¬ 
ing  the  cross-correlation  methpi  described  ^bogp.  The  line 
shows  the  chest  expansion  Aheasuf^d  by|t}i^pf%imiftipr  (n* 
belt,  scaled  and  shifted  to  nmUteh^fre^uficjyh^s^Thd  f'“' 
good  agreement  between  tlUse  dfeta^S^  inDifc^bs  maf fee 
frequency  shifts  were  most  likely  caused  by  respiration  I 
rather  than  by  cardiac  motion.  This  also  suggests  that  ^ 
physiologic  monitoring  of  the  chest  expansion  could  be  g 
used  to  correct  frequency  shifts  in  the  spectra. 

Frequency  Correction  Examples 

The  frequency  correction  method  is  demonstrated  in  Fig. 

3,  which  shows  a  series  of  single-shot  (NEX  =  1),  fixed  TE 
(TE  =  45  ms)  spectra  acquired  from  a  malignant  lesion. 

The  frequency  variation  evident  in  the  uncorrected  spectra 
(Fig.  3a)  is  removed  after  frequency  correction  (Fig.  3b). 

After  averaging  and  fitting  both  the  corrected  and  uncor¬ 
rected  series,  the  corrected  spectrum  (Fig.  4a,  top)  had 
improved  spectral  resolution,  increased  tCho  amplitude, 
and  decreased  fitting  error  as  compared  to  the  uncorrected 
spectrum  (Fig.  4a,  bottom).  This  example  is  not  represen¬ 
tative  of  a  typical  acquisition,  because  the  SNR  was  atyp- 
ically  high  and  TE  averaging  was  not  used.  However,  this 


FIG.  2.  Demonstration  of  the  relationship  between  spectral  fre¬ 
quency  shifts  and  chest  expansion.  The  subject  was  asked  to 
breathe  normally  while  a  series  of  localized,  unsuppressed  spectra 
were  acquired  with  TR  =  3  s.  The  line  shows  the  (smoothed)  chest 
expansion  as  measured  with  a  pneumatic  chest  bellows.  The  dots 
show  the  magnitude  of  frequency  shifts  measured  using  the  cross- 
correlation  method  on  the  spectra  at  corresponding  time  points.  The 
magnitude  of  the  chest  expansion  data  were  manually  scaled  to 
match  the  measured  frequency  offsets. 
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FIG.  3.  Demonstration  of  frequency  correc¬ 
tion.  A  series  of  64  identical  single-shot 
(NEX  =  1),  water-suppressed,  in  vivo  spec¬ 
tra  (TR/TE  -  3000/45  msec)  acquired  during 
normal  respiration  are  shown  (a)  prior  to  and 
(b)  after  frequency  correction  using  the 
cross-correlation  method.  The  measured 
frequency  shifts  In  this  acquisition  varied 
over  a  range  of  fl(Av)  with  a  SD  of  cr(Av). 
Note  that  these  data  are  not  TE-averaged 
and  have  an  unusually  high  signal-to-noise 
ratio. 


quency  correction.  Figures  4b  and  4c  show  how  the  mag¬ 
nitude  of  the  frequency  variation  can  affect  the  spectral 
quality.  Both  are  good-quality  spectra  with  resonance  from 
several  metabolites  in  addition  to  tCho.  In  Fig.  4b,  the 
corrected  and  uncorrected  spectra  appear  similar  because 
the  frequency  variation  was  relatively  small.  In  Fig.  4c, 
where  the  frequency  variation  was  larger,  the  corrected 
data  clearly  show  improved  spectral  resolution.  Fig.  4  days 
shows  a  case  where  tCho  was  detectable  with  frequency 
correction  and  not  detectable  (using  the  atCho  criterion) 
when  the  frequency  correction  was  omitted. 


Quantitative  Spectroscopy 

The  quantification  procedure  was  performed  using  fre¬ 
quency  correction  on  the  full  set  of  280  TE-a\^rage<iispec- 


a 


ppm 


tra  with  a  detectable  tCho  resonance.  These  spectra  were 
then  reprocessed  with  the  frequency  correction  step  re¬ 
moved.  The  resultant  spectral  broadening  generally  in¬ 
creased  the  Cram6r-Rao  errors,  making  the  tCho  resonance 
undetectable  in  69/280  (25%)  spectra  based  on  the  crtCho  < 
100%  criterion.  Of  the  remaining  211  spectra,  the  calcu¬ 
lated  [tCho]  values,  shown  in  Fig.  5a,  did  not  change  (P  = 
0.28,  two-tailed  t  test).  The  errors,  shown  in  Fig.  5b,  were 
~an  average  of  28%  larger  when  the  frequency  correction 
was  omitted.  As  expected,  the  increase  in  error  was  greater 
when  the  magnitude  of  the  frequency  variation  was  larger, 
as  shown  in  Fig.  5c.  The  range  and  SD  of  the  frequency 
Vfriation  are  shown  as  histograms  in  Fig.  5  days.  The  mean 

®ge  of  the  frequency  variation  R{ Av)  was  24  Hz 
4  ppm).  g, 

**£:■■■•.  fii 


FIG.  4.  Impact  of  frequency  correction  on  spectral  quantification.  Each  pair  of  spectra  shown  in  a-d  are  averaged  from  the  same  acquired 
data  but  processed  with  (top)  and  without  (bottom)  frequency  correction.  The  Inset  tables  show  calculated  [tCho]  concentrations  and  fitting 
errors  (mmol/kg)  based  on  the  method  In  Ref  14  for  both  corrected  an  uncorrected  spectra.  The  range  R( Av)  and  SD  o-(Av)  of  the  frequency 
variation  over  each  acquisition  is  also  shown.  The  pair  in  a  are  averaged  from  the  single-TE  arrays  shown  in  Figs.  3(a)  and  (b).  Pairs  b-d 
are  examples  of  TE-averaged  spectra  (TE  =  45-1 96  ms  in  1 28  increments,  TR  =  3  s),  demonstrating  cases  where  the  effect  of  frequency 
correction  is  relatively  small  (b),  moderate  (c),  and  large  (d).  In  the  uncorrected  spectrum  of  d,  the  normalized  fitting  error  (crtCho  =  1 40%) 
exceeded  the  detection  threshold,  so  in  this  case  a  tCho  resonance  Is  only  detectable  with  frequency  correction. 


FIG.  5.  Quantitative  impact  of 
frequency  correction  on  quantifi¬ 
cation  error  evaluated  in  a  series 
of  280  in  vivo  spectra,  (a)  Plot 
compares  the  total  tCho  concen¬ 
tration  ([tCho])  when  processed 
with  and  without  frequency  cor¬ 
rection.  There  is  no  significant  dif¬ 
ference  (P  =  0.28).  (b)  Plot  shows 
the  fitting  errors  (normalized  SD 
based  on  Cram6r-Rao  lower 
bounds)  when  processed  with 
and  without  frequency  correc¬ 
tion.  On  average,  the  error  Is  28% 
larger  without  frequency  correc¬ 
tion.  (c)  The  increase  in  error  is 
plotted  against  the  SD  of  Av  for 
each  spectrum,  (d)  Histograms 
for  the  range  R{ Av)  and  SD  <x(Av) 
of  the  frequency  variation  In  all 
280  spectra. 


DISCUSSION  f  |v  X 

This  study  demonstrates  that  the  correction  ot  respiratory- 
induced  frequency  variations  improves  the1  quality  of 
breast  MRS  at  4  T.  The  proposed  method  of  retrospec¬ 
tively  correcting  frequency  shifts  prior  to  averaging  is  Ap¬ 
plicable  in  cases  where  gross  displacement  of  the  tissue  is 
relatively  small,  but  B0  variations  are  substantial.  This 
scenario  is  typical  of  high-field  breast  MRS,  in  which 
breast  motion  is  restricted  by  the  coil  platform,  but  fre¬ 
quency  variations  are  large  dye  to  the  proximity  of  the 
lungs.  These  frequency  varidpBns  will  be  f^atSr.Bt  higb^. 
magnetic  fields  because  su&ep|ibi^ty-|ndSceK  fifidipisp 
tortions  scale  with  field  stjShgt®T%P^pfet^f  tMs  fe# 
fact  has  not  been  reported  at  1.5  T,  presumably  because  the 
susceptibility  effects  are  less  apparent  than  at  4  T. 

The  magnitude  of  the  respiration-induced  frequency 
variation  measured  in  the  breast  at  4  T  typically  ranges 
from  10  to  30  Hz,  with  a  mean  of  24  Hz  found  in  this  study. 
The  magnitude  depends  on  the  position  of  the  voxel 
within  the  breast.  All  the  Av  maps  acquired  in  this  study 
showed  a  general  trend  of  increased  frequency  variation 
closer  to  the  chest  wall.  Further  studies  are  required  to 
measure  the  full  spatial  dependence  of  Av  and  determine  if 
cardiac  motion  contributes  to  the  effect.  Other  factors  not 
explicitly  measured  but  expected  to  affect  the  artifact  in¬ 
clude  depth  of  respiration,  chest  size,  and  body  composi¬ 
tion.  If  left  uncorrected,  frequency  variations  during  an 
acquisition  will  cause  a  blurring  of  the  averaged  spectrum, 
reducing  the  effective  spectral  resolution  and  distorting 
the  line  shape. 

Although  the  qualitative  impact  of  frequency  shifts  on 
breast  spectra  is  often  not  dramatic,  particularly  when  the 


frequency  sHifts!are  small,  the  quantitative  impact  of  this 
artifact  is  abwaysjpresent.  Even  small  frequency  variations 
produce  afr  increase  in  fitting  errors  (Fig.  5c),  effectively 
decreasiril  the  sensitivity  of  spectroscopic  measurements, 
ifiis  iSMost  critical  when  fitting  metabolites  at  low  signal- 
tq-tioise  ratios,  #hich  is  common  in  breast  MRS. 

The  increased  error  is  likely  explained  by  the  increased 
dine  widths  produced  by  uncorrected  frequency  variation. 
In  general,  peak  amplitude  estimation  errors  are  larger  for 
broader  resonances.  This  has  beenshown  analytically  for 
Jhg-  casllfly^imple  j^prentffllpii  rgspiance  (24).  The  Voigt 
fiie-sha^OTodil  Jx|ed^iri|this  $ork  is  equivalent  to  a 
llorent^dn  linfl  shu^^cfey^Pe®  with  a  Gaussian  fre¬ 
quency  blurring,  which  enables  the  model  to  account  for 
frequency  variations  without  introducing  bias  in  the  am¬ 
plitude  estimates.  These  properties  were  verified  by  sim¬ 
ulating  spectral  data  sets  with  random  frequency  varia¬ 
tions:  as  the  magnitude  of  frequency  variation  increased, 
the  amplitude  estimates  did  not  change,  but  the  amplitude 
error  estimates  did  increase  (data  not  shown).  These  re¬ 
sults  are  consistent  with  the  in  vivo  data  shown  in  Fig. 
5a- c.  Although  not  evaluated  in  this  work,  it  is  expected 
that  other  spectral  fitting  methods  would  produce  similar 
results,  provided  a  sufficiently  flexible  line-shape  model  is 
used.  This  is  supported  by  a  previous  study  (11)  that 
showed  an  increase  in  fitting  errors  but  no  affect  on  peak 
amplitudes  when  fitting  spectra  with  the  LCModel  soft¬ 
ware  package  (25),  which  also  uses  a  flexible  line-shape 
model. 

There  are  a  number  of  approaches  for  handling  respira¬ 
tion-induced  frequency  variations.  In  principle,  the  arti¬ 
fact  can  be  avoided  altogether  using  respiratory  triggering 
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or  breath-hold  acquisitions.  These  methods  are  error- 
prone  and  are  less  efficient  in  use  of  machine  time,  but  are 
necessary  for  spectroscopic  editing  techniques  requiring 
precise  frequency  selectivity.  Alternatively,  the  frequency 
distortion  function  Av(f)  can  be  estimated  and  retrospec¬ 
tively  corrected  prior  to  averaging  spectra,  as  is  commonly 
done  for  phase  variations.  The  Av(f)  function  can  be  esti¬ 
mated  from  either  physiologic  monitoring,  MR  navigator 
signals,  or  directly  from  the  acquired  data.  Physiologic 
monitoring,  typically  done  with  chest  expansion  bellows 
or  ECGs,  is  least  direct  and  requires  careful  processing  to 
properly  denoise  and  detrend  the  data.  Using  navigator 
signals  has  the  disadvantage  that  there  is  typically  a  delay 
between  the  navigator  acquisition  and  the  actual  FID;  the 
length  of  the  delay  may  be  significant  when  compared  to 
the  respiratory  period.  Although  not  shown,  we  also  de¬ 
veloped  and  tested  a  version  of  LASER  with  an  interleaved 
small-flip-angle  STEAM  acquisition  immediately  follow¬ 
ing  the  metabolite  acquisition  period,  similar  to  the 
method  described  by  Thiel  et  al.  (7).  This  also  worked 
acceptably,  but  it  was  not  as  effective  as  the  cross-correla¬ 
tion  method,  probably  due  to  the  delay  between  the  navi¬ 
gator  echo  and  the  measured  FID  (350  ms).  A  nonlocalized 
navigator  with  a  shorter  delay  was  also  tested,  but  these 
shift  measurements  did  not  correlate  well  with  the  local¬ 
ized  shifts.  f 

Extracting  Av(f)  directly  from  the  data,  as  in  this' work, Is 
the  most  direct  way  of  measuring  the  artifact;  This  ap¬ 
proach  requires  sufficient  SNR  in  each  acquisition  to  mea¬ 
sure  frequency  shifts.  Numerous  methods  ha^e  been  pro¬ 
posed  for  measuring  Av  between  consecutive  spectra,  in¬ 
cluding  cross-correlation,  frequency-referencing  a  single  . 
resonance,  time  domain  fitting,  and  principal  comporfi§ntv 
analysis  (18,26).  The  time  domain  method  does  not  work 
well  with  TE-averaged  data  because  the  initial  portion  of 
each  individual  FID  is  distorted  by  B0  modulations,  whiefi  ‘ 
are  corrected  only  after  averaging.  The  single  resonance 
referencing  method  (with  the  1.3-ppm  lipid)  and  cross¬ 
correlation  method  worked  c||bparably,  bufithp,  cross-cor¬ 
relation  method  was  overall Hore  fthtfit^pfeausei^t^S 
the  full  spectral  informatioaf.y^jsi^ .  jpe  vgal 

resonance  as  a  frequency  reference  fi  as  b  e en  prop  os e d  by 
other  groups  (10,12).  This  approach  is  also  feasible  pro¬ 
vided  the  water  suppression  is  adjusted  properly  and  has 
sufficient  bandwidth  to  avoid  frequency-dependent  line- 
shape  variations  like  those  shown  in  Fig.  3a.  Methods 
using  principle  component  analysis  may  also  be  effective, 
but  their  applicability  for  TE-averaged  spectra  has  not  yet 
been  investigated. 

This  study  did  not  measure  spectroscopic  phase  varia¬ 
tions  due  to  physiologic  motion.  Although  in  studies  of 
compliant  subjects  there  were  generally  no  large-scale  mo¬ 
tions  of  the  breast,  small  displacements  of  the  breast  (par¬ 
ticularly  near  the  chest  wall)  were  visible  by  imaging. 
These  small-scale  breast  motions  may  produce  phase  vari¬ 
ations  (10).  With  the  TE-averaged  acquisition  method  used 
in  this  work,  the  zero-order  phase  varies  with  each  TE  due 
to  B0  modulations,  so  phase  correction  is  required  even  in 
the  absence  of  motion.  Therefore,  separating  physiologic 
and  system  sources  of  phase  variation  was  not  feasible. 


CONCLUSION 

This  study  reports  measurements  of  respiratory-induced 
frequency  variations  in  the  breast  at  4T  and  their  impact 
on  quantitative  MRS.  The  average  frequency  variation 
measured  over  a  typical  acquisition  was  24  Hz,  which  is 
approximately  10X  greater  than  comparable  measure¬ 
ments  of  this  effect  in  the  brain.  If  left  uncorrected,  this 
artifact  reduces  the  ability  to  measure  metabolite  concen¬ 
trations.  In  our  sample  of  280  spectra  with  a  detectable 
tCho  resonance,  quantitative  fitting  errors  increased  by  an 
average  of  28%,  and  in  69  (25%)  spectra  the  tCho  reso¬ 
nance  became  undetectable.  These  results  indicate  that 
using  frequency  correction  can  improve  the  quality  of 
quantitative  breast  1H  MRS  at  high  field. 
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In  Vivo  Quantification  of  Choline  Compounds  in  the 
Breast  With  1H  MR  Spectroscopy 

Patrick  J.  Bolan,1,2  Sina  Meisamy,1,2  Eva  H.  Baker,1,2  Joseph  Lin,1,2  Timothy  Emory,2 
Michael  Nelson,2  Lenore  I.  Everson,2  Douglas  Yee,3,4  and  Michael  Garwood1-3* 


This  work  describes  a  methodology  for  quantifying  levels  of 
total  choline-containing  compounds  (tCho)  in  the  breast  using 
in  vivo  1H  MR  spectroscopy  (MRS)  at  high  field  (4  Tesla).  Water 
is  used  as  an  internal  reference  compound  to  account  for  the 
partial  volume  of  adipose  tissue.  Peak  amplitudes  are  esti¬ 
mated  by  fitting  one  peak  at  a  time  over  a  narrow  frequency 
band  to  allow  measurement  of  small  metabolite  resonances  in 
spectra  with  large  lipid  peaks.  This  quantitative  method  signif¬ 
icantly  improves  previously  reported  analysis  methods  by  ac¬ 
counting  for  the  variable  sensitivity  of  breast  1H  MRS  measure¬ 
ments.  Using  this  technique,  we  detected  and  quantified  a  tCho 
peak  in  214  of  500  in  vivo  spectra.  tCho  levels  were  found  to  be 
significantly  higher  in  malignancies  than  in  benign  abnormali¬ 
ties  and  normal  breast  tissues,  which  suggests  that  this  tech¬ 
nique  could  be  used  to  diagnose  suspicious  lesions  and  mon¬ 
itor  response  to  cancer  treatments.  Magn  Reson  Med  50: 
1134-1143,  2003.  ©  2003  Wiley-Liss,  Inc. 

Key  words:  choline;  breast  cancer;  quantification;  MRS;  diag¬ 
nosis 

Breast  cancer  is  a  very  common  disease,  affecting  11%  of 
American  women  and  causing  more  than  40000  deaths 
each  year  (1).  While  breast  cancer  mortality  is  decreasing, 
the  incidence  continues  to  rise  (2).  Thus,  there  is  a  great 
need  for  noninvasive  diagnostic  tools  for  both  screening 
and  treatment  monitoring.  The  conventional  diagnostics— 
X-ray  mammography,  sonography,  and  physical  examina¬ 
tion — are  limited  in  their  sensitivity  for  detecting  disease 
and  their  specificity  for  distinguishing  between  benign 
and  malignant  lesions.  Magnetic  resonance  imaging  (MRI) 
of  the  breast  is  being  used  increasingly  because  of  its  high 
sensitivity,  but  its  reported  specificity  is  widely  variable  (3). 

Researchers  have  recently  begun  to  augment  breast  MRI 
studies  with  MR  spectroscopy  (MRS)  to  increase  specific- 
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ity.  In  vivo  MRS  can  detect  a  resonance  at  3.25  ppm  that 
has  contributions  from  several  different  compounds,  in¬ 
cluding  choline,  phosphocholine,  glycerophosphocholine, 
and  taurine.  High-resolution  in  vitro  and  ex  vivo  studies 
indicate  that  the  levels  of  choline  compounds  increase 
with  malignancy  (4-6).  At  the  lower  field  strengths  used 
for  in  vivo  work  (1.5-4  T),  these  multiple  resonances 
cannot  be  spectrally  resolved  and  thus  appear  as  a  single 
peak,  termed  total  choline-containing  compounds  (tCho). 

Several  studies  conducted  at  1.5  T  have  shown  that  in 
vivo  MRS  can  be  used  to  distinguish  between  benign  and 
malignant  tissues  (7-11).  These  studies  used  the  hypoth¬ 
esis  that  tCho  is  only  detectable  in  malignancies.  A  pooled 
analysis  of  these  five  studies  showed  that  this  tCho  detect¬ 
ability  criterion  can  identify  malignancies  with  an  83% 
sensitivity  and  85%  specificity  (12).  This  qualitative  ap¬ 
proach  is  promising,  but  it  is  only  applicable  if  the  MRS 
measurement  sensitivity  is  invariant.  In  similar  studies 
performed  at  4  T,  the  increased  sensitivity  allows  detec¬ 
tion  of  tCho  in  benign  lesions  and  normal  subjects.  A  more 
general  approach  is  to  quantify  the  tCho  peak  with  the 
expectation  that  tCho  levels  are  higher  in  malignancies 
than  in  benign  lesions  or  normal  tissues.  Two  groups  have 
reported  quantification  of  tCho  levels  using  external  phan¬ 
tom  referencing  methods  (7,13).  These  studies  demon¬ 
strated  the  feasibility  of  quantitative  breast  MRS,  but  they 
were  limited  to  small  patient  groups  and  did  not  report 
measurement  errors.  Quantification  is  also  valuable  for 
measuring  tumor  response  to  treatment  regimens.  Two 
other  studies  found  that  the  detectability  and  amplitude 
(judged  qualitatively)  of  the  tCho  peak  decreases  after  che¬ 
motherapy  (8,11).  A  quantitative  method  enables  finer 
measurements  of  the  magnitude  and  rate  of  tumor  re¬ 
sponse. 

Although  quantification  of  metabolite  levels  is  routinely 
performed  in  MRS  of  the  brain,  it  is  more  difficult  to 
perform  in  the  breast  because  of  the  heterogeneous  distri¬ 
bution  of  the  glandular  and  adipose  tissues.  Although  the 
spectroscopist  will  typically  plan  a  voxel  to  include 
mainly  glandular  tissue  and  tumor,  voxels  of  typical  size 
(1-2  mL)  nearly  always  contain  some  adipose  tissue  as 
well.  The  amount  of  included  adipose  tissue  can  vary 
greatly  depending  on  the  architecture  of  the  gland  and/or 
lesion. 

The  two  basic  elements  of  a  quantitative  MRS  method¬ 
ology  are  the  referencing  strategy  and  the  spectral  fitting 
technique.  The  referencing  strategy  proposed  in  this  work 
uses  water  as  an  internal  reference  peak.  This  approach 
compensates  for  the  partial  volume  of  adipose  tissue  in  the 
voxel  and  naturally  leads  to  a  molal  (mol/kg)  concentra¬ 
tion  for  water-soluble  metabolites.  The  fitting  technique 
used  in  this  work  is  based  on  a  hybrid  time-domain  (TD) 
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and  frequency-  domain  (FD)  method  called  “TDFD  fitting” 
(14).  This  method  enables  flexible  lineshape  definition  by 
the  use  of  a  TD  model,  and  has  excellent  frequency-selec¬ 
tion  properties  since  the  residuals  are  evaluated  and  min¬ 
imized  in  the  FD,  The  ability  to  select  a  narrow  frequency 
range  is  crucial  for  fitting  small  resonances  in  the  presence 
of  very  large  ones,  as  is  the  case  in  breast  spectra  contain¬ 
ing  large  lipid  peaks. 

The  goal  of  this  project  was  to  develop  a  method  to 
quantitatively  measure  tCho  levels  in  breast  tissue.  The 
methodology  presented  integrates  several  existing  tech¬ 
niques:  single-voxel  localization  by  adiabatic  selective  re¬ 
focusing  (LASER)  (15),  TE  averaging  to  reduce  lipid  side¬ 
band  artifacts  (16),  automatic  frequency  referencing  to  cor¬ 
rect  respiration  artifacts,  frequency-selective  spectral 
fitting,  and  quantification  using  water  as  an  internal  refer¬ 
ence  peak.  We  quantified  500  in  vivo  spectra  with  this 
technique,  analyzed  the  results,  and  compared  them  with 
an  independent  method  based  on  external  referencing. 
The  applicability  of  using  quantitative  1H  MRS  for  diag¬ 
nosing  suspect  breast  lesions  is  discussed  below. 

MATERIALS  AND  METHODS 

Acquisition 

All  measurements  were  performed  with  a  hybrid  4  T  sys¬ 
tem,  consisting  of  a  90-cm-bore  magnet  (model  4  T-900; 
Oxford  Magnet  Technology,  Oxfordshire,  UK)  with  a  clin¬ 
ical  gradient  system  (Sonata;  Siemens,  Erlangen,  Germany) 
interfaced  with  an  imaging  spectrometer  (Unity  Inova; 
Varian,  Palo  Alto,  CA).  Several  different  single-breast 
quadrature  transmit/receive  RF  surface  coils  of  similar 
design  were  used  to  accommodate  different  breast  sizes. 
The  coils  were  mounted  on  a  custom-built  patient  table 
designed  for  unilateral,  prone  breast  studies. 

A  total  of  105  subjects  (23-72  years  old,  mean  = 
48  years)  were  studied  in  175  MRI/MRS  sessions.  Of  these, 
86  were  participants  in  a  study  involving  the  diagnosis  of 
lesions  with  suspicious  mammographic  findings,  14  were 
participants  in  a  study  regarding  monitoring  response  to 
neoadjuvant  chemotherapy,  and  five  were  presumed  nor¬ 
mal  volunteers  (no  breast-related  health  problems  or  ab¬ 
normal  mammograms).  Approximately  half  of  these  stud¬ 
ies  were  performed  after  a  needle  biopsy  or  other  invasive 
procedure.  All  studies  were  approved  by  the  institutional 
review  board,  and  informed  written  consent  was  obtained 
from  the  subjects  prior  to  the  studies. 

All  100  patients  from  the  diagnosis  and  treatment-mon¬ 
itoring  studies  were  examined  with  the  same  MRI/MRS 
protocol,  which  consisted  of  high-resolution  imaging,  dy¬ 
namic  contrast-enhanced  imaging,  and  single-voxel  spec¬ 
troscopy.  The  subjects  were  positioned  prone  with  their 
breast  centered  horizontally  in  the  magnet.  After  scout 
images  were  acquired  to  verify  position,  the  coil  was  man¬ 
ually  tuned  and  matched.  A  high-resolution  3D  fast  low- 
angle  shot  (FLASH)  image  (fat-suppressed,  matrix  -  256  x 
256  X  64,  field  of  view  (FOV)  =  14-18  cm,  TE/TR  = 
4.1/13.5  ms,  flip  angle  =  30°)  and  a  fast  2D  multislice 
image  (fat-suppressed,  matrix  =  256  X  128,  30  slices,  slice 
thickness  =  2.5  mm,  FOV  =  14-18  cm,  TE/TR  =  5.1/390  ms, 
flip  angle  =  90°)  were  acquired  prior  to  injection  of  Gd-DTPA 


(0.1  mmol/kg  body  weight).  Five  to  seven  2D  images  were 
acquired  immediately  after  injection,  followed  by  a  second 
3D  image.  Both  image  sets  were  analyzed  using  own  image- 
processing  software  (developed  with  Matlab  (The  Math- 
works,  Natick,  MA))  to  select  voxels  for  MRS  with  the  subject 
still  in  the  magnet.  The  criteria  for  voxel  selection  included 
lesion  architecture,  dynamic  Gd-DTPA  uptake,  and  prior 
clinical  information  from  mammographic  or  ultrasound  im¬ 
ages.  The  voxels  were  planned  to  maximize  coverage  of  the 
lesion  and  minimize  the  inclusion  of  adipose  tissue.  The  five 
normal  subjects  were  studied  with  only  the  high-resolution 
3D  FLASH  image  followed  by  spectroscopy,  since  no  Gd- 
DTPA  was  administered  in  these  studies. 

All  spectra  were  acquired  using  the  LASER  localization 
technique  (15)  with  4096  complex  points  and  6-kHz  spec¬ 
tral  width.  Each  voxel  measurement  began  with  a  calibra¬ 
tion  of  the  localized  R,  field  strength,  followed  by  30-60  s 
of  manual  adjustment  of  the  linear  shims.  A  fully  relaxed, 
single-shot,  unsuppressed  spectrum  was  acquired  to  mea¬ 
sure  the  water  and  lipid  peaks.  The  power  required  to 
suppress  the  water  signal  using  variable  pulse  power  and 
optimized  relaxation  delays  (VAPOR)  (17)  water  suppres¬ 
sion  was  then  manually  adjusted.  The  metabolite  spec¬ 
trum  was  acquired  using  TE  averaging  with  TE  =  45-196 
ms  in  64  or  128  increments  and  with  TR  =  3  s  (16).  Each 
free  induction  decay  (FID)  signal  was  individually  saved, 
and  no  averaging  was  performed  until  processing.  Each 
voxel  required  —9  min  in  total.  One  to  four  voxels  were 
studied  in  each  subject,  for  a  total  acquisition  time  of 
~1  hr. 

Preprocessing  Spectra 

All  spectral  processing  programs,  including  preprocess¬ 
ing,  fitting,  and  quantification,  were  written  in  Matlab  (The 
Mathworks,  Natick,  MA).  In  all  spectra,  the  last  512  points 
of  the  raw  TD  FID  were  used  to  calculate  the  root-mean- 
square  (RMS)  noise  and  correct  for  DC  offsets.  For  the 
unsuppressed  spectra  used  to  measure  the  water  and  lipid 
peaks,  the  FIDs  were  truncated  to  1024  points  and  then 
zero-filled  to  2048  points.  The  zero-order  phase  4>  was 
measured  and  corrected  using  the  average  of  two  different 
autophasing  methods:  1)  fitting  the  phase  of  the  first  few 
TD  points  to  a  line  and  using  4>  at  time  t  ~  0,  and  2)  finding 
the  value  of  <|)  that  maximizes  the  smallest  value  of  the  real 
part  of  the  FD  spectrum.  Using  the  mean  of  these  two 
methods  produces  a  robust  estimate  of  4>.  The  spectra  were 
then  frequency  referenced  by  setting  the  maximum  of  the 
water  peak  to  4.7  ppm.  For  the  water-suppressed  spectra 
used  to  measure  the  metabolites,  each  individual  spectrum 
from  the  TE  averaging  acquisition  was  automatically 
phased  as  described  above.  To  correct  respiration-induced 
frequency  shifts,  each  spectrum  was  shifted  in  frequency 
to  maximize  the  cross-correlation  function  between  it  and 
the  first  spectrum  of  the  acquisition  (18).  After  phasing  and 
frequency  correction  was  performed,  the  spectra  were  av¬ 
eraged. 

Fitting 

To  fit  the  water,  1.3  ppm  lipid,  and  tCho  peaks  in  the 
processed  spectra,  a  new  fitting  program  was  developed. 
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suppressed  and  averaged  spectra,  tCho  was  fit  at  3.25  ppm. 
The  baseline  around  3.25  ppm  is  often  distorted,  so  a 
linear  baseline  model  was  used  to  reduce  bias.  Figure  1 
shows  the  relationship  between  the  data  (Sdata),  model 
(Smode!),  baseline  (Sbaselin0),  and  residual  (fl).  Using  a  lin- 
ear  baseline,  the  residual  function  to  be  minimized  is 

fl(«)  =  [sdaj<o)  -  -  ^Sn™w(0}K;S:5SS!  [2] 

where  u)0  is  the  center  of  the  frequency  range  being  fit.  No 
baseline  correction  was  used  for  fitting  the  water  and  lipid 
peaks,  so  in  those  cases  Sbasoline  =  0  for  all  w. 

The  peak  amplitude  A  is  the  parameter  of  greatest  inter¬ 
est  since  it  is  proportional  to  the  number  of  nuclei  in  the 
voxel.  The  fitting  error  was  estimated  using  the  Cramer- 
Rao  minimum  variance  bound  (CRMVB)  of  the  parameter 
A: 


FIG.  1.  Model  showing  the  frequency-selective  fit  of  a  tCho  peak 
with  baseline  correction.  The  real  component  of  the  raw  data  is  fit 
with  a  Voigt  lineshape  model  (thick  line)  plus  a  linear  baseline 
(dashed  line).  The  residual  (dotted  line)  is  minimized  over  a  0.4-ppm 
region  centered  around  the  peak.  The  same  model  without  the  linear 
baseline  is  used  to  fit  the  unsuppressed  water  and  lipid  peaks. 


The  algorithm  is  based  on  the  TDFD  fitting  technique 
proposed  by  Slotboom  et  al.  (14).  In  this  approach,  a  TD 
model  is  used  to  describe  the  peak  parameters,  but  the 
fitting  is  performed  by  applying  a  Fourier  transform  and 
minimizing  the  residuals  in  the  FD.  This  flexible  fitting 
method  has  excellent  frequency-selective  properties  and 
can  be  used  with  various  lineshape  definitions.  In  this 
implementation,  a  Voigt  lineshape  was  used  to  model  all 
peaks  (19).  The  modeled  signal  smodei  from  a  single  peak  s 
can  be  described  in  the  TD  t  as 

Smodei(f)  =  A  *  exp(-itof  +  i'4>  -  Xt  -  y2f)  [1] 

with  signal  amplitude  A ,  chemical  shift  frequency  o),  zero- 
order  phase  <(>,  Lorentzian  damping  factor  X,  and  Gaussian 
damping  factor  y.  The  Lorentzian  and  Gaussian  damping 
factors  are  related  to  the  full  width  at  half  maximum 
(FWHM)  of  an  impulse  response:  X  =  rr  •  FWHM  and  y  = 
tt  ■  FWHMI (2V\n  2),  as  described  by  Ogg  et  al.  (20).  Note 
that  although  all  equations  in  this  section  are  shown  with 
continuous  time  and  frequency  for  clarity,  all  calculations 
were  performed  with  discrete  variables. 

The  peaks  were  fit  one  at  a  time  by  minimizing  the  FD 
residuals  over  a  0.4-ppm  (68-Hz)  section  of  the  spectrum 
centered  on  the  peak.  The  zero-order  phase  <|>  was  fixed  at 
zero  (since  it  was  presumably  corrected  during  preprocess¬ 
ing)  and  only  the  real  portion  of  the  residual  was  mini¬ 
mized.  Initial  values  for  parameters  were  generated  based 
on  prior  knowledge  and  simple  heuristics  (e.g.,  the  initial 
frequency  is  that  of  the  maximum  of  the  absolute  value 
spectrum  over  the  fitting  region).  The  nonlinear  minimi¬ 
zation  of  the  residuals  was  performed  using  the  optimiza¬ 
tion  toolkit  provided  with  Matlab  (function  Isqnonlin ,  us¬ 
ing  the  large-scale  model  option). 

In  the  unsuppressed  spectra,  water  was  fit  at  4.7  ppm, 
and  the  polymethylene  lipid  peak  was  fit  at  1.3  ppm.  In  the 


CRMVBa  - - -  [31 

where  a£oiSB  is  the  variance  of  white,  Gaussian  noise.  The 
CRMVB A  is  not  actually  an  estimate  of  the  fitting  error; 
rather,  it  is  a  theoretical  minimum  limit  for  the  estimation 
accuracy.  Nevertheless,  the  CRMVB A  is  commonly  used  in 
NMR  applications  to  estimate  fitting  errors  (21).  Effec¬ 
tively,  this  measure  is  the  noise  variance  scaled  by  the 
sensitivity  of  the  model  to  changes  in  the  parameter  A.  The 
error  is  often  expressed  as  a  normalized  standard  deviation 
(SD):  crA  =  V CRMVB A/ A. 

This  error  estimate  was  also  used  to  establish  the  detec¬ 
tion  criteria.  In  all  cases,  a  fit  was  initially  performed  and 
the  parameter  estimates  and  errors  were  calculated.  If  the 
fitting  error  for  A  was  greater  than  a  specified  threshold, 
then  the  fit  was  rejected  and  the  resonance  was  considered 
undetectable.  For  all  of  the  spectra  reported  in  this  paper, 
the  threshold  for  the  normalized  error  was  unity:  if  crA  >  1, 
the  peak  was  considered  undetectable.  Although  this 
threshold  is  statistically  arbitrary,  it  is  convenient  and 
roughly  corresponds  to  a  signal-to-noise  ratio  (SNR)  of 
2-3. 


Quantification 

The  spectral  fitting  produces  an  amplitude  A  for  each 
peak,  expressed  in  arbitrary  units  (au).  To  standardize  this 
measurement  for  use  in  both  internal  and  external  refer¬ 
encing  schemes,  several  corrections  must  be  made  to  ac¬ 
count  for  experimental  conditions  and  the  physical  prop¬ 
erties  of  each  species.  The  corrected  amplitude  A '  is: 


/gain  =  gain/gain0 

[5] 

/coii  =  B,/S,p0 

[6] 
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/t,  «  1  -  exp(— TR/T,) 

17] 

.  N 

/tj  “  ^  S  exp(-TE/T2), 

[8] 

j-i 


where  gain  is  the  receiver  gain,  is  the  local  amplitude  of 
the  excitation  radiofrequency  field,  TR  is  the  pulse  repeti¬ 
tion  time,  and  TE  is  the  echo  time.  The  receiver  correction 
factor  /gain  is  necessary  if  the  receiver  gain  is  different  in 
the  suppressed  and  unsuppressed  acquisitions.  The  coil 
receive  efficiency  factor  /coil  was  calculated  by  assuming 
the  transmit  and  receive  efficiencies  are  equal.  The  refer¬ 
ence  values  gain0  and  are  arbitrary,  but  must  be  used 
consistently  when  comparing  values  from  different  acqui¬ 
sitions.  The  Tz  correction  /Tt  is  approximate,  and  is  valid 
provided  TR  >>  TE.  With  a  TE-averaged  acquisition,  the 
T2  correction/^  is  a  summation  over  N acquisitions,  each 
with  different  a  TE  (16).  Because  it  is  impractical  to  mea¬ 
sure  relaxation  properties  in  each  voxel,  constant  values 
were  assumed  for  all  relaxation  constants  based  on  mea¬ 
surements  in  several  subjects. 

After  these  corrections  are  made,  the  signal  amplitudes 
are  proportional  to  the  number  of  nuclei  in  the  volume. 
The  ratio  of  the  tCho  and  water  amplitudes  can  be  con¬ 
verted  to  molal  concentration  (moles  solute  per  mass  sol¬ 
vent)  by  correcting  for  the  number  of  1H  nuclei  per  mole¬ 
cule  n  and  the  molecular  weight  of  the  solvent  MWwater: 

[tcho] = t9] 

Note  that  the  coil  efficiency  factor  /coil  cancels  because 
both  water  and  tCho  come  from  the  same  volume  of  inter¬ 
est.  This  quantity  [tCho],  expressed  in  units  mmol/kg,  is 
the  metric  proposed  as  an  in  vivo  measure  of  the  tissue 
level  of  choline-containing  compounds  in  the  breast.  This 
measurement  is  presented  along  with  the  SD  of  the  fitting 
error:  [tCho]  ±  [tCho]  •  <xA. 

In  spectra  where  no  tCho  peak  was  detected  based  on  the 
crA  >  1  criterion,  an  additional  procedure  was  performed 
to  determine  the  sensitivity  of  the  measurement.  A  simu¬ 
lated,  noiseless  spectrum  containing  a  single  peak  at 
3.25  ppm  with  a  Gaussian  FWHM  of  15  Hz  and  amplitude 
Asfm  was  added  to  the  preprocessed  in  vivo  spectrum.  The 
combined  in  vivo  +  simulated  spectrum  was  then  fit  using 
the  procedure  described  above,  and  the  fitting  error  <rA  was 
calculated.  This  process  was  repeated  with  successively 
smaller  simulated  peak  amplitudes,  reduced  10%  each 
iteration,  until  the  simulated  peak  was  no  longer  detect¬ 
able.  The  smallest  value  of  Asim  that  led  to  a  detectable 
tCho  peak  was  then  used  to  calculate  a  minimum  detect¬ 
able  level  (MDL)  of  tCho,  [tCho]MDL,  using  Eq.  [9J.  For 
these  in  vivo  spectra  where  tCho  is  not  detectable,  the 
overall  [tCho]  measurement  is  expressed  as  0  ±  [tCho]MDL. 

To  validate  this  proposed  internal  referencing  scheme, 
the  tCho  level  was  also  calculated  independently  using  an 
external  referencing  scheme.  In  general,  the  corrected  sig¬ 
nal  amplitude  A'  of  a  resonance  is  proportional  to  the 
number  of  nuclei  n  in  the  sample:  n  =  Ksy3A'.  The  system 


constant  Ksys  (with  units  mol/au)  accounts  for  the  system- 
specific  hardware  and  software.  The  value  of  Ksys  was 
calculated  for  our  system  in  separate  calibration  experi¬ 
ments  with  a  phantom  of  known  concentration  and  mea¬ 
surable  relaxation  properties.  The  externally-referenced 
concentration  of  tCho  can  then  be  expressed  in  molal  units 
(mol /kg)  as 

[tCh°L'  =  (ftMJr}lCh0  ^  [101 

where  pwator  is  the  water  density,  and  V  is  the  voxel  vol¬ 
ume.  To  compare  the  internal  and  external  methods,  the 
tissue  water  density  was  assumed  to  be  1  kg/L.  The  volume 
V  was  assumed  to  be  the  entire  voxel  volume,  ignoring  the 
effect  of  partial  volume  from  adipose  tissue. 

RESULTS 

The  quantification  method  was  applied  to  500  spectra 
acquired  from  105  subjects  in  175  MR  study  sessions.  The 
voxel  size  ranged  from  0.4  to  16  mL,  with  a  median  volume 
of  1.6  mL.  A  peak  near  3.25  ppm  was  detected  in 
214  (43%)  of  the  500  spectra.  The  processing  and  fitting 
processes  required  ~1  min  per  spectrum  on  a  typical 
workstation.  After  we  ran  the  fully-automated  procedure 
on  all  500  spectra  and  manually  reviewed  the  results,  six 
spectra  required  reprocessing  with  manual  intervention  to 
adjust  phase  and/or  frequency  referencing.  Each  water, 
lipid,  and  tCho  peak  was  fit  with  a  Voigt  lineshape,  giving 
both  Lorentzian  and  Gaussian  contributions  to  the  line- 
width.  The  fit  water  peak  had  a  median  Lorentzian  and 
Gaussian  linewidth  of  12  Hz  and  13  Hz,  respectively;  the 
1.3  ppm  lipid  peak  had  16  Hz  and  28  Hz,  and  the  tCho 
peak  had  0  Hz  and  14  Hz. 

To  calculate  absolute  tCho  levels,  relaxation  values  were 
measured  in  several  subjects  using  standard  spectroscopic 
techniques.  The  measured  values  with  SDs  were: 
^i.wator  “  870  ~  325  ms,  T2  wator  —  60  ±  7  ms,  — 

480  ±  100  ms,  T2  lipid  =  69  ±  12  ms,  and  T2  tCho  =  399  ± 
133  ms.  Tltcho  was  not  measured,  but  was  assumed  to  be 
the  same  as  Tl  Avater.  Although  some  of  these  SDs  are  large, 
their  overall  effect  on  the  [tCho]  measurement  is  relatively 
small.  Using  our  acquisition  parameters,  the  uncertainties 
in  T-i. water,  T2.watsr,  T,  tch0,  and  T2,tCho  led  to  [tCho]  SDs  of 
0%,  5%,  6%,  and  13%,  respectively,  for  a  total  of  16%. 

Several  examples  of  in  vivo  spectra  are  shown  in  Fig.  2. 
Figure  2a  shows  a  large  voxel  acquired  in  healthy  glandu¬ 
lar  tissue.  A  clear  tCho  peak  is  visible  at  3.25  ppm,  and  the 
fitting  produces  a  measurement  of  [tCho]  =  0.66  ± 
0.06  mmol/kg.  The  model  fit  of  this  peak  is  shown  above 
the  full  spectrum,  and  the  smooth  residual  is  shown  be¬ 
neath.  There  is  another  metabolite  peak  visible  at  3.4  ppm 
that  was  not  fit.  This  sample  shows  that  tcho  can  be 
detected  in  normal  breast  tissue.  Figure  2b  is  a  spectrum 
from  an  invasive  ductal  carcinoma  with  ItCho]  =  6.1  ± 
0.08  mmol/kg.  The  volume  of  this  voxel  is  smaller  than 
that  of  the  voxel  in  Fig.  2a,  but  since  it  is  located  closer  to 
the  coil  the  sensitivity  is  comparable.  The  SNR  of  the  tCho 
peak  is  quite  high,  but  the  residual  shows  some  structure, 
indicating  an  imperfect  fit.  Figure  2c  shows  a  more  typi- 
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FIG.  2.  Example  spectra.  A  water-suppressed  spectrum  is  shown  on  the  right,  with  the  tCho  fit  shown  above  and  the  residual  (including 
linear  baseline)  underneath.  The  location  of  the  voxel  is  shown  in  a  contrast-enhanced,  fat-suppressed  sagittal  image  on  the  left,  a:  Normal 
gland,  volume  =  13.0  mL,  (tCho]  =  0.66  ±  0.06  mmol/kg,  lipid  fraction  =  3%.  b:  Malignant  tumor  of  invasive  ductal  carcinoma,  volume  = 
6.8  mL,  [tCho]  =  6.1  ±  0.08  mmol/kg,  lipid  fraction  -  8%.  c:  Benign  finding  of  atypical  hyperplasia  In  an  Insensitive  region  of  the  coil, 
volume  =1.1  mL,  [tCho]  =  1.4  ±  0.7  mmol/kg,  lipid  fraction  =  14%. 


cally-sized  voxel  in  a  lesion  later  identified  by  needle 
biopsy  as  atypical  hyperplasia,  with  [tCho]  =  1.5  ± 
0.8  mmol/kg.  Atypical  hyperplasia  is  generally  considered 
to  be  benign,  but  it  is  a  marker  indicating  an  increased  risk 
for  developing  future  malignancies.  The  low  SNR  of  the 
tCho  peak  is  reflected  in  the  error  estimated  by  the  fitting 
procedure. 

These  examples  demonstrate  how  the  sensitivity  of 
breast  MRS  can  vary  greatly.  Due  to  variability  in  coil 
loading,  voxel  size,  and  partial  volume  of  adipose  tissue, 
the  sensitivity  for  detecting  tCho  varied  by  a  factor  of  >100 
in  this  study.  Figure  3  shows  the  fitting  error  in  all 
500  spectra  as  a  function  of  water  SNR  and  voxel  volume. 


For  both  natural  and  simulated  spectra,  the  fitting  error  is 
the  normalized  error  multiplied  by  the  calculated  concen¬ 
tration,  [tCho]  *  <rA  or  [tCho]MDL  *  <rA.  Clearly,  the  water 
SNR  is  a  better  indicator  of  fitting  error,  since  it  automat¬ 
ically  corrects  for  coil  efficiency  and  the  partial  volume  of 
adipose  tissue.  These  plots  show  reasonable  properties  for 
an  unbiased  fitting  method.  For  example,  fitting  error  de¬ 
creases  uniformly  with  increasing  water  SNR,  and  tCho  is 
more  likely  to  be  detected  in  spectra  with  lower  fitting 
errors. 

Although  the  fitting  error  is  greater  in  smaller  voxels,  the 
[tCho]  measurement  itself  is  independent  of  the  voxel  size, 
as  shown  in  Fig.  4a.  The  filled  diamonds  represent  [tCho] 


Quantification  of  tCho  in  Breast 


1139 


0.1  1  10  100 
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FIG.  3.  Factors  contributing  to  the  sensitivity  of  [tCho]  measure¬ 
ment.  a:  The  fitting  error  from  all  500  spectra  as  a  function  of  the 
SNR  of  the  unsuppressed  water  peak  fTD  amplitude/RMS  noise),  b: 
The  fitting  error  as  a  function  of  total  voxel  volume.  Riled  circles  are 
error  estimates  from  spectra  where  tCho  was  detectable;  open 
circles  are  error  estimates  calculated  using  the  MDL  procedure  (see 
text)  from  spectra  where  no  tCho  was  detectable. 


measurements  in  spectra  where  a  peak  was  detectable,  and 
the  hollow  diamonds  represent  the  MDL  of  tCho  in  spectra 
where  no  peak  was  detectable.  As  expected,  the  smaller 
voxels  were  less  likely  to  have  detectable  tCho,  but  clearly 
the  voxel  size  did  not  bias  the  [tCho]  measurement  in 
spectra  with  detectable  tCho. 

Figure  4b  shows  how  the  [tCho]  measurement  varies 
with  the  lipid  content  of  the  voxel.  The  lipid  fraction  is 
estimated  using  the  ratio  between  the  corrected  ampli¬ 
tudes  of  the  water  and  1.3  ppm  lipid  peaks:  %  -  A;ipid/ 
Mwutor  +  A  ijpid)  X  100.  The  [tCho]  measurement  was 
expected  to  be  independent  of  the  voxel  lipid  fraction 
because  all  of  the  metabolites  that  contribute  to  the  tCho 
peak  are  water-soluble.  In  voxels  with  low  to  moderate 
lipid  content,  [tCho]  was  independent  of  the  lipid  fraction. 
In  voxels  with  large  lipid  content,  however,  the  [tCho] 


measurement  increases  with  increasing  lipid  fraction.  Lip¬ 
ids  apparently  contribute  to  the  amplitude  of  the  tCho 
peak,  either  through  baseline  artifacts  not  suppressed  by 
TE  averaging  or  by  a  true  resonance  at  3.25  ppm.  Because 
of  this  lipid  contamination,  we  arbitrarily  chose  a  cutoff 
value  of  33%,  above  which  the  [tCho]  measurement  is 
considered  biased  and  the  MRS  measurement  is  consid¬ 
ered  invalid. 

Figure  5  shows  a  comparison  between  the  internal  and 
external  reference  schemes  described  in  Materials  and 
Methods.  This  shows  that  in  good-quality  spectra,  where 
tCho  is  detectable  and  the  lipid  content  is  low,  both  inter¬ 
nal  and  external  methods  produce  consistent  results  (R2  - 
0.91).  The  absolute  values  produced  by  the  external 
method  are  lower,  due  to  overestimation  of  the  product  (V* 
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FIG.  4.  Dependence  of  [tCho]  on  voxel  volume  and  lipid  content. 
The  [tCho]  measurements  from  spectra  where  tCho  was  detectable 
are  Indicated  by  solid  diamonds;  for  spectra  where  tCho  did  not 
meet  the  detectability  threshold,  the  MDLs  of  tCho  are  indicated  by 
hollow  diamonds.  Part  a  shows  that  the  [tCho]  measurement  is 
independent  of  voxel  volume,  which  is  desirable.  Part  b  shows  that 
when  the  lipid  volume  fraction  is  low,  the  [tCho]  measurement  is 
unbiased.  When  the  lipid  content  Is  high,  the  [tCho]  measurement 
increases.  Spectra  with  lipid  fractions  >33%  (an  arbitrarily  selected 
threshold)  are  considered  biased. 
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FIG.  5.  A  comparison  between  [tCho]  calculated  with  both  internal 
and  external  referencing  methods  in  98  spectra,  all  of  which  had 
detectable  tCho  and  a  lipid  fraction  of  ^33%.  A  linear  fit  through  the 
origin  gives  R2  =  0.91  and  a  slope  of  0.72.  This  shows  that  these 
measures  are  highly  correlated,  although  the  external  method  pro¬ 
duces  somewhat  smaller  values. 


p yvatBr)  in  the  denominator  of  Eq.  [10],  which  represents  the 
aqueous  content  of  the  voxel. 

The  applicability  of  the  [tCho]  measurement  for  diagnos¬ 
ing  different  pathologies  is  shown  in  Fig.  6.  All  spectra 
that  were  of  sufficient  quality  (i.e.,  free  of  artifacts,  lipid 
fraction  <  33%)  and  were  acquired  from  lesions  with 
biopsy-confirmed  pathology  were  divided  into  a  "malig¬ 
nant”  category  (including  infiltrative  ductal,  lobular,  and 
unspecified  carcinomas)  and  a  “benign”  category  (includ¬ 
ing  atypical  hyperplasias,  fibroadenomas,  fibrocystic 
changes,  and  cysts).  There  were  insufficient  data  to  distin¬ 
guish  further  histological  subcategories  or  tumor-staging 
grades.  When  multiple  spectra  were  acquired  from  a  single 
lesion,  only  the  spectrum  with  the  smallest  error  was 
included  in  this  chart.  Spectra  labeled  "normal”  were  ac¬ 
quired  from  normal  volunteers  and  from  regions  of  nor¬ 
mal-appearing  and  asymptomatic  glandular  tissue  in  other 
subjects.  These  results  show  that  the  tCho  measurement  is 
elevated  in  malignancies  and  some  benign  lesions.  The 
mean  [tCho]  is  greater  in  malignancy  than  in  benign  tis¬ 
sues  (P  =  0.008,  one-tailed  f-test),  but  the  difference  be¬ 
tween  the  normal  and  benign  categories  is  not  statistically 
significant  (P  =  0.17).  An  ROC  analysis  was  performed  to 
determine  a  threshold  [tCho]  value  for  distinguishing  be¬ 
tween  benign  and  malignant  lesions.  Using  equal  weight¬ 
ing  for  false  positives  and  false  negatives,  the  criteria  for 
malignancy  is  [tCho]  >  1.38  mmol/kg.  With  this  cutoff,  the 
sensitivity  is  46%  and  the  specificity  is  94%.  Note  that 
neither  of  these  analyses  take  into  account  the  variable 
sensitivity. 

Single-voxel  MRS  in  the  breast  is  very  sensitive  to  the 
size  and  placement  of  the  voxel  because  of  the  heteroge¬ 
neous  distribution  of  tCho  in  the  breast,  as  can  be  seen  in 
the  spectra  shown  in  Fig.  7.  All  three  spectra  were  ac¬ 
quired  from  different  regions  of  the  same  3-cm  tumor,  a 
grade  III  invasive  ductal  carcinoma  that  was  studied  after 


the  patient  had  received  4  months  of  chemotherapy.  In  a 
large  voxel  covering  most  of  the  tumor,  [tCho]  was  mea¬ 
sured  to  be  1.1  ±  0.6  mmol/kg.  The  anterior,  enhancing 
portion  of  the  tumor  had  a  higher  [tCho]  (1.5  ± 
0.5  mmol/kg),  whereas  the  posterior,  non-enhancing  re¬ 
gion  had  no  detectable  peak  ([tCho]  —  0  ±  1.3  mmol/kg). 
The  sensitivity  is  lowest  in  the  posterior  voxel  due  to  its 
distance  from  the  coil,  small  voxel  size,  and  higher  lipid 
content.  This  example  underscores  the  importance  of 
proper  voxel  placement. 

DISCUSSION 

This  work  describes  a  new  method  for  quantifying  tCho 
levels  in  breast  tissue.  Quantitative  MRS  is  a  substantial 
improvement  over  the  qualitative  detection  methods  used 
in  previous  studies  of  breast  MRS.  Quantification  is  par¬ 
ticularly  important  in  the  breast  because  the  sensitivity  of 
the  MRS  measurement  is  generally  more  variable  than  it  is 
in  brain  tissue.  This  is  due  primarily  to  the  highly  variable 
adipose  tissue  content  of  the  breast  and  the  greater  varia¬ 
tion  of  the  coil  receive  efficiency.  Previous  in  vivo  studies 
used  the  hypothesis  that  detectability  of  tCho  is  associated 
with  malignancy;  however,  this  approach  is  only  valid  if 
the  detection  threshold  is  constant.  We  have  found  that  the 
use  of  a  higher  B0  field  (4  T)  and  optimized  surface  coils 
increases  the  sensitivity  enough  to  enable  the  detection  of 
tCho  in  normal  breast  tissue  and  several  benign  lesions. 

Two  previous  reports  used  external  referencing  to  quan¬ 
tify  choline  levels  in  breast  tissue  (7,13).  Neither  of  these 
studies  corrected  for  partial  volume  of  adipose  tissue 
within  the  voxel.  In  our  experience,  the  amount  of  intra¬ 
voxel  adipose  tissue  can  vary  greatly,  as  shown  in  Fig.  4b. 
The  breast  can  be  coarsely  approximated  by  a  two-com- 
partment  model  consisting  of  aqueous  regions  (fibroglan- 
dular  tissue)  and  aliphatic  regions  (adipose  tissue).  The 
aqueous  compartment  contains  all  of  the  choline-contain¬ 
ing  compounds  that  are  known  to  be  elevated  in  malignan- 


FIG.  6.  Measurements  of  [tCho]  In  malignant,  benign,  and  normal 
tissues.  Error  bars  represent  SDs  in  spectra  where  tCho  was  de¬ 
tectable,  and  MDLs  in  spectra  where  no  tCho  was  detectable. 
Measurements  are  ordered  from  the  largest  to  the  smallest  in  each 
category.  Malignant  and  benign  spectra  were  from  biopsy-con- 
firmed  lesions  only.  Normal  spectra  were  selected  from  healthy- 
appearing  glandular  tissue. 
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FIG.  7.  Spatial  variation  of  [tCho]  measurements  within  a  single  tumor.  All  spectra  are  from  a  3-cm  tumor  of  invasive  ductal  carcinoma  after 
4  months  of  chemotherapy  treatment.  Images  on  the  left  show  voxel  positions  on  a  sagittal  slice  from  a  contrast-enhanced,  fat-suppressed 
3D  FLASH  Image.  Water-suppressed  spectra  on  the  right  are  shown  with  the  tCho  fit  above  and  the  residual  (including  baseline) 
underneath.  In  spectrum  C,  where  no  tCho  was  detectable,  the  MDL  of  tCho  is  shown  with  a  dotted  line,  a:  Voxel  containing  both  enhancing 
and  non-enhancing  regions  of  the  tumor,  volume  =  7.8  mL,  [tCho)  =  1 .1  ±  0.6  mmol/kg,  lipid  fraction  =  5%.  b:  Voxel  acquired  from  the 
most  enhancing  region  of  the  tumor,  volume  =  1.1  mL,  [tCho]  -  1.5  ±  0.5  mmol/kg,  lipid  fraction  =  2%.  c:  Voxel  acquired  from  a 
non-enhancing  region,  volume  =  0.9  mL,  [tCho]  =  0  ±  1.3  mmol/kg,  lipid  fraction  -  13%. 


cies.  Any  externally -referenced  method  must  account  for 
this  compartmentation  to  avoid  systematic  error,  whereas 
an  internally-referenced  method  compensates  for  compart¬ 
mentation  automatically.  Internal  referencing  methods 
have  further  advantages  in  that  they  do  not  require  sepa¬ 
rate  calibration  experiments,  and  they  use  fewer  error- 
introducing  correction  factors  for  coil  loading,  coil  effi¬ 
ciency,  and  voxel  volume. 

The  disadvantage  of  the  internal  referencing  approach  is 
that  it  requires  the  density  of  NMR-visible  water  in  the 
aqueous  compartment  to  be  relatively  constant.  There  are 
physiological  conditions  that  affect  the  water  content,  but 
these  effects  are  expected  to  be  small  compared  to  the 


variability  of  compartmentation.  Both  internal  and  exter¬ 
nal  approaches  will  be  affected  by  the  presence  of  edema 
from  previous  invasive  procedures  or  radiation  therapy, 
and  by  variations  in  cellularity  and  relaxation  rates.  Pre¬ 
vious  reports  expressed  concern  that  the  relaxivity  of  gad¬ 
olinium  contrast  agents  would  deleteriously  affect  MRS 
measurements  (8).  In  our  experience,  gadolinium  effects 
are  small  ( — 10%)  compared  to  the  measurement  errors 
and  variability  of  tCho  levels  (22). 

The  spectral  fitting  software  used  in  this  work  is  highly 
specialized  for  fitting  peaks  in  breast  spectra,  and  is  avail¬ 
able  from  the  authors  by  request.  The  method  is  based  on 
the  TDFD  fitting  technique,  which  uses  an  analytical  TD 
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model  but  minimizes  residuals  in  the  FD  (14).  We  modi¬ 
fied  that  technique  by  adding  a  linear  baseline  correction 
and  using  a  very  narrow  frequency  band  (0.4  ppm)  to  fit  a 
single  peak  at  a  time.  By  using  only  the  real  component  of 
the  spectrum,  along  with  the  narrow  frequency  band  and 
the  linear  baseline  correction,  small  tCho  peaks  can  be  fit 
without  bias  from  large  neighboring  lipid  resonances.  An¬ 
other  useful  feature  of  the  TDFD  method  is  that  it  permits 
the  use  of  the  Voigt  lineshape,  which  can  fit  purely  Lorent- 
zian  or  purely  Gaussian  lineshapes,  or  any  combination  of 
the  two.  This  flexibility  was  useful  for  this  application 
because  the  use  of  a  Lorentzian  or  Gaussian  model  alone 
did  not  produce  good-quality  fits  for  all  peaks.  The  tCho 
peak  is  difficult  to  model  precisely  because  it  is  a  super¬ 
position  of  several  resonances,  and  at  4  T  there  is  not 
sufficient  spectral  resolution  to  separate  the  individual 
components.  Also,  imperfect  correction  of  respiration-in¬ 
duced  frequency  shifts  can  cause  “blurring"  that  leads  to 
Gaussian  lineshapes  in  averaged  spectra.  Although  the 
Voigt  model  did  not  always  produce  an  ideal  fit  (e.g.,  Fig. 
2b),  it  performed  reasonably  well  in  most  cases. 

The  TDFD  method  was  chosen  because  it  is  well  suited 
to  the  problem  and  is  relatively  easy  to  implement.  There 
are  numerous  other  methods  that  could  be  adapted  for 
fitting  breast  spectra;  however,  most  of  these  methods  are 
incapable  of  fitting  a  precise  frequency  range  or  analyti¬ 
cally  modeling  a  non-Lorentzian  lineshape  (23,24).  The 
LCModel  method  (25)  is  very  useful  for  fitting  brain  spec¬ 
tra,  and  it  may  also  be  applicable  for  fitting  breast  spectra. 
With  the  currently  available  implementation  of  the  LC¬ 
Model  software,  we  were  unable  to  get  good-quality  fits  of 
small  tCho  peaks  in  spectra  with  large  lipids.  If  the  LC¬ 
Model  method  could  be  adapted  for  breast  spectra,  it  might 
offer  some  advantages,  such  as  producing  more  accurate 
fits  due  to  its  sophisticated  baseline  model,  and  enabling 
the  fitting  of  other  metabolites  that  are  occasionally  ob¬ 
servable  in  breast  spectra,  such  as  creatine,  glycine,  lac¬ 
tose,  and  taurine.  Automated  fitting  methods  are  generally 
desirable  because  they  eliminate  user  interaction  and  pro¬ 
vide  a  well  accepted  measure  of  error  (Cramer-Rao 
bounds)  unbiased  by  peak  lineshape.  However,  a  simple 
analysis  using  peak  integration  and  linear  baseline  correc¬ 
tion  over  manually-defined  extents  gave  comparable  re¬ 
sults  (data  not  shown). 

Investigators  commonly  use  CRMVBs  to  estimate  the 
variance  of  model  parameters  when  fitting  NMR  spectra 
(21).  In  this  work,  the  fitting  error  is  used  as  an  estimate  of 
the  overall  measurement  error,  although  it  is  possibly  an 
underestimate.  Most  notably,  the  CRMVB  errors  from  fit¬ 
ting  the  water  are  very  small,  but  the  actual  error  in  esti¬ 
mating  the  water  amplitude  is  higher.  This  occurs  because 
the  Cramer-Rao  theory  assumes  that  the  data  are  perfectly 
described  by  the  model  function,  and  this  assumption  is 
violated  by  irregular  water  lineshapes  and  the  overlapping 
5.4  ppm  lipid  resonance.  Additionally,  certain  experimen¬ 
tal  factors  (e.g.,  patient  motion  and  respiration)  can  make 
the  true  measurement  error  greater  than  the  fitting  error. 
The  true  measurement  error  can  only  be  established 
through  repeatability  studies.  Nevertheless,  some  estimate 
of  the  measurement  error  is  essential  for  any  quantification 
scheme,  and  it  should  be  incorporated  into  the  interpreta¬ 
tion  of  quantitative  MRS  results.  It  is  particularly  impor- 
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tant  in  breast  MRS  because  of  its  highly  variable  sensitiv¬ 
ity,  as  demonstrated  in  Fig.  3a.  Note  that  this  sensitivity 
can  be  estimated  prior  to  the  acquisition  of  an  entire  spec¬ 
trum  by  using  the  SNR  of  the  water  peak  in  a  single-shot, 
unsuppressed  spectrum.  This  information  can  be  used  to 
prescreen  voxels  for  a  diagnostic  study. 

The  iterative  MDL  method  described  here  (in  which  a 
simulated  tCho  peak  is  added  to  a  spectrum  and  fit)  can  be 
used  to  measure  the  sensitivity  in  spectra  with  no  detect¬ 
able  tCho  peak.  This  is  critical  for  interpreting  whether  a 
negative  finding  of  tCho  is  due  to  lack  of  sensitivity  or  low 
tCho  levels.  The  MDL  is  closely  related  to  the  fitting  error: 
under  ideal  conditions,  the  fitting  error  is  equal  to  the 
MDL,  but  in  practice  the  MDL  is  often  larger. 

The  use  of  a  quantitative  method  to  measure  tCho  levels 
in  the  breast  increases  the  usefulness  of  MRS  for  diagnos¬ 
ing  benign  and  malignant  lesions.  Figure  6  shows  [tCho] 
measurements  divided  into  three  broad  pathological  cate¬ 
gories.  The  finding  that  tCho  levels  are  higher  in  malig¬ 
nancies  than  in  benign  or  normal  tissues  is  consistent  with 
previously  published  reports  that  used  the  detection  of 
tCho  to  indicate  malignancy  (7-11).  The  ItCho]  threshold 
of  1.38  mmol/kg,  based  on  an  ROC  analysis,  can  be  used  to 
distinguish  between  malignant  and  benign  lesions;  how¬ 
ever,  the  sensitivity  is  poor  (46%).  More  data  must  be 
acquired  to  determine  whether  subdividing  the  malignant 
and  benign  categories  will  improve  these  results  (e.g.,  in 
the  benign  category,  the  two  highest  values  were  atypical 
hyperplasias).  Note  also  that  this  ROC  analysis  is  only  a 
first-order  method  of  interpreting  the  [tCho]  levels.  More 
sophisticated  analyses  should  be  used  to  account  for  the 
measurement  errors  and  [tChol  probability  distribution 
functions. 

The  absolute  tCho  levels  (0.4-10  mmol/kg)  reported  in 
this  work  are  reasonably  consistent  with  previous  in  vivo 
estimates.  Roebuck  et  al.  (7)  found  choline  levels  of  0.4- 
5.8  mmol/L,  and  estimated  their  detection  threshold  was 
0.2  mmol/L.  Bakken  et  al.  (13)  reported  a  single  measure¬ 
ment  of  2  mmol/L.  Further  investigations  are  required  to 
determine  whether  the  absolute  [tCho]  values  reported 
here  are  reproducible  on  different  MR  scanners  and  in 
other  institutions. 

Figure  6  shows  a  large  range  of  [tChol  measurements 
in  known  cancers,  from  0.5  to  >8  mmol/kg.  tCho  levels 
did  not  appear  to  be  related  to  different  histological 
types  of  cancer  (e.g.,  lobular  vs.  ductal);  however,  there 
were  insufficient  data  to  establish  statistical  signifi¬ 
cance.  The  variation  in  [tCho]  was  also  not  explained  by 
complicating  factors  such  as  recent  biopsy  or  previous 
chemotherapy  treatments,  although  some  of  the  spectra 
in  which  no  tCho  was  detected  were  acquired  from 
tumors  that  had  been  recently  biopsied.  The  large  range 
in  [tCho]  measurements  may  be  a  natural  feature  of 
cancer,  due  to  variations  in  the  amount  of  water,  lipids, 
and  stroma  within  the  lesion  and  adjacent  tissues,  and 
to  variations  in  intracellular  tCho  levels  and  neoplastic 
cell  density.  Previous  chemotherapy,  radiation,  or  inva¬ 
sive  procedures  may  further  complicate  these  issues. 
Additionally,  certain  experimental  factors  can  contrib¬ 
ute  to  inaccuracy  in  these  measurements — most  notably 
patient  motion  and  incorrect  voxel  positioning.  If  either 
the  biopsy  needle  or  the  MRS  voxel  is  placed  in  the 
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wrong  region,  there  will  be  no  agreement  between  the 
[tCho]  measurements  and  the  histology.  Finally,  it 
should  be  noted  that  MRS  measurements  can  be  ad¬ 
versely  affected  by  metallic  radiographic  markers, 
which  are  widely  used  but  cannot  always  be  identified 
in  heterogeneous  breast  tissue. 

CONCLUSIONS 

In  this  work  we  report  a  new  technique  for  quantifying 
tCho  levels  in  breast  cancer  using  single-voxel  1H  MRS. 
With  the  use  of  optimized  surface  coils  and  a  high-field 
(4  Tesla)  scanner,  we  measured  tCho  levels  in  normal 
breast  tissue  and  in  benign  and  malignant  lesions.  The 
levels  of  tCho  were  found  to  be  elevated  in  malignancies 
compared  to  benign  lesions,  indicating  that  quantitative 
MRS  may  be  used  to  aid  in  diagnosing  breast  lesions  and 
monitoring  response  to  cancer  treatments. 
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Abstract 


Purpose:  To  determine  if  changes  in  the  concentration  of  choline-containing  compounds 
([tCho])  from  before  to  within  24  hours  after  receiving  primary  systemic  therapy  (PST) 
predict  clinical  response  in  patients  with  locally  advanced  breast  cancer  (LABC). 

Materials  and  Methods:  Sixteen  women  with  biopsy-confirmed  LABC  scheduled  for 
doxorubicin-based  PST  were  recruited.  Four  Tesla  MRI/MRS  scans  were  performed  prior 
to  treatment,  within  24  hours  after  the  first  dose,  and  after  the  fourth  dose  of  PST. 

Lesions  size  was  assessed  with  Gd— DTPA  enhanced  MRI.  Lesion  [tCho]  was  quantified 
using  single-voxel  proton  ('H)  MRS.  Statistical  analyses  used  Pearson’s  correlation 
coefficient  and  Wilcoxon  rank-sum  test. 

Results:  Fourteen  of  16  patients  completed  the  protocol.  [tCho]  was  not  measurable  in 
one  patient  due  to  unfavorable  lesion  morphology  for  MRS  voxel  placement.  Of  the 
remaining  13  patients,  4  had  inflammatory  breast  cancer,  6  had  invasive  ductal  carcinoma 
(IDC),  2  had  invasive  lobular  carcinoma  (ILC),  and  1  had  mixed  IDC/ILC.  Based  on 
RECIST  criteria,  8/13  patients  were  objective  responders  and  5/13  were  nonresponders. 
The  change  in  [tCho]  between  baseline  and  within  24  hours  after  the  first  dose  of  PST 
showed  significant  positive  correlation  with  the  change  in  lesion  size  (R  =  0.79,  p  = 
0.001).  The  change  in  [tCho]  within  24  hours  was  significantly  different  between 
objective  responders  and  nonresponders  (p  =  0.007). 


2 


Conclusion:  These  results  suggest  that  the  change  in  [tCho]  between  the  baseline  scan 
and  within  24  hours  after  the  first  dose  of  PST  can  serve  as  an  indicator  for  predicting 
clinical  response  to  doxorubicin-based  chemotherapy  in  LABC. 


Introduction 


Primary  systemic  therapy  (PST),  also  known  as  neoadjuvant  chemotherapy,  given 
prior  to  breast  cancer  surgery  offers  several  advantages  over  standard  post-operative 
chemotherapy.  Although  PST  does  not  offer  any  survival  benefits  over  post-operative 
chemotherapy  (1,2),  patients  who  receive  PST  are  more  likely  to  receive  breast- 
conserving  surgery.  Moreover,  the  use  of  PST  permits  in  vivo  monitoring  of  tumor 
response.  Complete  disappearance  of  tumor  at  surgical  resection  (pathologic  complete 
response)  was  associated  with  the  best  overall  survival  in  multiple  studies  (2-4).  Since 
there  are  many  active  agents  to  treat  breast  cancer,  it  is  important  to  know  early  in  the 
course  of  treatment  whether  the  drug  chosen  will  be  effective  for  an  individual. 

Currently  there  are  no  standardized  criteria  that  can  individually  detect  early 
response  to  PST.  Conventional  modalities  such  as  physical  examination,  ultrasonography, 
and  mammography  vary  in  reliability  for  measuring  tumor  response  (5-8).  MRI  is 
increasingly  being  used  to  evaluate  locally  advanced  breast  cancer  (LABC)  defined  as 
invasive  carcinomas  >  2  cm  in  longest  diameter  with  or  without  inflammatory  features. 
With  respect  to  treatment  monitoring,  studies  have  shown  correlation  between  specific 
MRI  findings  and  clinical  response.  Changes  in  lesion  size,  neoplastic  phenotype, 
dynamic  contrast  enhancement,  and  extraction  flow  product  all  correlate  with  clinical ' 
response  (9-15).  However,  changes  in  lesion  size  or  dynamic  contrast  enhancement 
measured  by  MRI  are  not  detected  until  at  least  6  weeks  following  PST  (13).  The  ability 
to  immediately  detect  response  to  a  specific  chemotherapeutic  regimen  would  be  ideal, 
since  it  would  allow  for  optimal  individualization  of  chemotherapeutic  regimens  for 
patients  with  the  goal  of  obtaining  a  pathologic  complete  response. 
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Recently  there  has  been  an  interest  in  the  use  of  proton  magnetic  resonance 
spectroscopy  (*H  MRS)  for  detecting  and  monitoring  patients  with  breast  cancer.  Using 
‘H  MRS,  it  has  been  shown  that  neoplastic  breast  tissue  contains  elevated  levels  of 
choline-containing  compounds  (tCho)  which  yield  a  signal  at  a  chemical  shift  of  3.2  ppm 
(16-20).  One  group  of  researchers  has  shown  that  changes  in  the  tCho  signal  occur  in 
breast  cancer  patients  who  receive  PST  (21).  Other  studies  have  shown  a  correlation 
between  the  change  in  tCho  signal  and  clinical  response  in  patients  with  extracranial 
lymphomas  and  germ  cell  tumors  (22)  and  pediatric  gliomas  (23). 

'ft  MRS  of  the  breast  is  technically  challenging  because  sensitivity  can  be 
limiting  and  spectral  artifacts  can  occur.  Due  to  the  heterogeneous  distribution  of  fat  and 
glandular  tissue  in  the  breast,  *H  spectra  of  breast  often  contain  large  lipid  signals  that 
give  rise  to  contaminant  peaks  around  3.2  ppm.  Fortunately,  these  artifactual  peaks,  also 
known  as  sidebands,  can  be  suppressed  using  a  recently  developed  technique  called  echo¬ 
time  (TE)  averaging  (24).  Additionally,  treatment  monitoring  requires  implementation  of 
a  method  to  quantify  spectra.  Previous  studies  at  lower  magnetic  field  strength  (1.5  T) 
were  based  on  the  hypothesis  that  tCho  is  detectable  only  in  malignant  breast  tissue  (20). 
However,  with  the  increased  sensitivity  at  high  field  (4T),  it  has  been  shown  that  tCho  is 
also  detectable  in  benign  lesions  and  normal  fibroglandular  breast  tissue  (25). 

We  hypothesized  that  an  early  decrease  in  the  concentration  of  choline-containing 
compounds  ([tCho])  could  identify  an  early  response  to  PST.  Thus,  the  purpose  of  our 
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study  was  to  determine  if  changes  in  [tCho]  from  before  to  within  24  hours  after  PST 
predict  clinical  response  in  patients  with  LABC. 

Materials  and  Methods 
Patients 

Women  between  the  ages  of  18  —  80  years  old  with  biopsy-confirmed  LABC  who 
were  scheduled  to  receive  doxorubicin-based  PST  were  eligible  to  enroll  in  our 
prospective  study.  Patients  were  enrolled  in  our  study  from  August  2001  to  April  2003. 
The  study  protocol  was  approved  by  the  Institutional  Review  Board  at  the  University  of 
Minnesota,  School  of  Medicine.  Informed  written  consent  was  obtained  prior  to  all 
studies.  Patients  were  referred  by  medical  oncologists  working  at  the  University  of 
Minnesota  or  in  private  practice.  Information  collected  for  each  patient  included  age, 
menopausal  status,  use  of  hormone  replacement  therapy,  use  of  oral  contraceptive  pills, 
method  of  breast  biopsy,  diagnosis,  ER  (estrogen  receptor)/PR  (progesterone  receptor) 
status,  lesion  size,  and  [tCho],  Each  patient  received  the  standard  4  doses  of  adriamycin 
and  cyclophosphamide  (AC)  (doxorubicin  —  60  mg/m2,  cyclophosphamide  —  600 
mg/m2).  The  first  dose  of  AC  was  given  on  day  1  and  each  additional  dose  of  AC  was 
given  at  2 1  -day  intervals  for  a  total  of  64  days.  Patients  underwent  MRI/MRS  scans  prior 
to  treatment,  within  24  hours  after  the  first  dose  of  AC  was  administered  (day  2),  and 
after  the  fourth  dose  of  AC  prior  to  surgery.  Patients  had  the  option  of  returning  for 
additional  scans  after  the  second  and  third  doses  of  AC.  Patients  received  a  baseline  scan 
at  a  median  of  2  days  (range  1  to  21  days)  prior  to  the  first  dose  of  AC  and  the  final  scan 
at  a  median  of  3  days  (range  1  to  8  days)  after  the  fourth  dose  of  AC.  After  completing  4 
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cycles  of  AC,  one  patient  continued  treatment  with  4  cycles  of  paclitaxel  and  underwent 
additional  scans  after  the  second  and  fourth  dose  of  paclitaxel. 

MRI  and  MRS  techniques 

All  measurements  were  performed  with  a  4T  research  scanner,  consisting  of  a 
90cm  bore  magnet  (model  4T  —  900,  Oxford  Magnet  Technology,  Oxford  shire,  UK) 
with  a  clinical  gradient  system  (model  'Sonata',  Siemens,  Erlangen,  Germany)  interfaced 
with  an  imaging  spectrometer  (model  'Unity  Inova',  Varian,  Palo  Alto,  CA).  The  gradient 
system  was  capable  of  40mT/m  with  a  400ps  rise  time.  Several  different  single-breast 
quadrature  transmit/receive  radio  frequency  (RF)  surface  coils  of  similar  design  were 
used  to  accommodate  different  breast  sizes  (26).  The  coils  were  mounted  onto  a  custom- 
built  patient  table  designed  for  unilateral,  prone  breast  studies.  Patients  were  positioned 
prone  with  their  breast  centered  horizontally  in  the  magnet.  After  acquiring  scout  images 
to  verify  proper  positioning,  the  coil  was  manually  tuned  and  matched.  A  high  spatial 
resolution  3D  fast  low-angle  shot  (FLASH)  image  (fat-suppressed,  matrix  =  256  x  256  x 
64,  field  of  view  =  14—18  cm,  TR/TE  =  13.5/4.1  ms,  flip  angle  =  30°)  and  a  fast  2D 
multislice  FLASH  image  (fat-suppressed,  matrix  =  256  x  128, 30  slices,  slice  thickness  = 
2.5  mm,  field  of  view  =  14—18  cm,  TR/TE  =  390/5.1  ms,  flip  angle  =  90°)  were  acquired 
prior  to  injection  of  Gd— DTP  A  (0.1  mmol/kg  body  weight).  Five  sets  of  2D  images  were 
acquired  immediately  after  injection,  followed  by  a  second  3D  image  acquisition.  Both 
image  sets  were  analyzed  using  our  own  image  processing  software  (developed  with 
Matlab;  The  Mathworks,  Natick,  MA)  to  select  voxels  for  MRS  with  the  subject  still  in 
the  magnet. 
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Voxel  placement  was  usually  performed  jointly  by  at  least  two  of  the  authors 
(SM,  PJB,  EHB,  and  MG).  In  a  few  instances  only  the  physician  (SM  or  EHB)  performed 
voxel  placement,  but  only  after  having  at  least  6  months  of  experience  in  breast  MR. 
Criteria  for  voxel  selection  included  lesion  architecture,  lesion  size  (>  2  cm  in  longest 
diameter),  dynamic  contrast  uptake,  and  prior  clinical  information  obtained  from 
mammographic  or  ultrasound  images.  Voxels  were  planned  to  maximize  coverage  of  the 
lesion  while  minimizing  inclusion  of  adipose  tissue.  Voxels  that  were  planned  in  studies 
following  PST  were  placed  to  cover  the  same  portion  of  the  lesion  covered  in  the  baseline 
scan.  If  the  lesion  changed  in  size,  the  voxel  was  adjusted  accordingly.  Single-voxel  'H 
spectroscopy  was  performed  with  the  technique  known  as  localization  by  adiabatic 
selective  refocusing  (LASER)  (27).  MRS  data  were  acquired  using  4096  complex  points 
and  6  kHz  spectral  width.  Each  voxel  measurement  began  with  a  calibration  of  the 
localized  Bi  field  strength,  followed  by  30—60  s  of  manual  adjustment  of  the  linear 
shims.  A  fully  relaxed,  single-shot,  unsuppressed  spectrum  was  acquired  to  measure  the 
water  and  lipid  signals.  The  RF  power  was  then  manually  adjusted  to  suppress  the  water 
signal  using  the  variable  pulse  power  and  optimized  relaxation  delays  (VAPOR) 
technique  (28).  The  metabolite  spectrum  was  acquired  using  TE  averaging  with  TE  = 
45—196  ms  in  64  or  128  increments  with  TR=3  s  (24).  Each  free  induction  decay  signal 
was  individually  saved  -  no  averaging  was  performed  until  processing.  Levels  of  [tCho] 
were  quantified  by  fitting  a  Voigt  lineshape  model  to  the  data  and  using  the  unsuppressed 
water  signal  as  an  internal  reference  (25).  All  [tCho]  measurements  were  recorded  in 
units  of  mmol/kg  water. 

Data  and  Statistical  Analysis 
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Statistical  Analysis  was  done  using  the  SAS  software  package  system  for 
Windows  (version  8.02,  Cary,  North  Carolina).  Lesion  size  was  measured  from  the  high 
spatial  resolution  3D  subtraction  image  that  was  obtained  by  subtracting  the  baseline  3D 
image  from  the  contrast-enhanced  3D  image.  The  variable  used  to  characterize  lesion  size 
was  longest  diameter  (LD).  LD  measured  at  baseline  and  after  the  fourth  dose  of  AC  are 
represented  by  the  variables  LDo  and  LDf  respectively.  Derived  variables  included  LD 
after  the  fourth  dose  of  AC  normalized  to  baseline  and  expressed  as  a  percentage, 


%LDf  =100 


[l] 


and  the  percent  change  in  lesion  size  from  baseline  to  after  the  fourth  dose  of  AC, 

%ALDf  =  lOof  (LDf_~LDo)l  r 

V  J  1 

Concentrations  of  tCho  at  baseline,  within  24  hours  after  the  first  dose  of  AC,  and  after 
the  fourth  dose  of  AC  are  represented  by  the  variables  [tCho]o,  [tCho]24,  and  [tCho]f, 
respectively.  Derived  variables  included  the  [tCho]  at  24  hours  normalized  to  baseline 
and  expressed  as  a  percentage, 

o/  _  i  c\r\(  [tCh°]  24 

/o[tCho]24  - 100^  [tCho]o  J  ,  [ 

and  the  percent  change  in  [tCho]  from  baseline  to  24  hours  after  the  first  dose  of  AC, 

•/„A[tCho]24=ioo(»i^2Lj  , 
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Using  the  classification  system  of  the  Response  Evaluation  Criteria  in  Solid 
Tumors  (RECIST)  (29),  patients  were  categorized  into  one  of  two  groups:  objective 
responders  or  nonresponders.  Objective  responders  were  patients  who  had  %ALDf  <  - 
30%.  Nonresponders  were  patients  who  had  %ALDf  >  -30%.  The  association  between 
%A[tCho]24  and  %ALDf  was  calculated  using  Pearson's  correlation  coefficient  (30).  The 
comparison  between  %A[tCho]24  and  patient  response,  and  the  comparison  of  [tChojo 
between  objective  responders  and  nonresponders,  were  calculated  using  the  Wilcoxon 
rank-sum  test  (3 1). 

Results 

Patient  and  Tumor  Characteristics 

A  total  of  16  patients  were  recruited,  of  which  14  completed  the  protocol.  Two 
patients  did  not  complete  the  study  because  one  opted  to  have  surgery  prior  to  completing 
all  four  doses  of  PST  and  the  second  patient  decided  to  undergo  an  alternative  method  of 
treatment.  Among  the  14  patients  who  completed  the  study,  1  patient  was  not  included  in 
the  data  analysis  due  to  an  inability  to  plan  a  suitable  voxel  for  [tCho]  measurements. 

This  patient  was  a  49  year-old  female  with  inflammatory  breast  cancer  (IBC).  At 
baseline,  the  maximum  width  of  the  mass  was  only  0.3  cm,  whereas  LDo  was  6  cm. 

Table  1  lists  the  lesion  size,  [tCho],  diagnosis,  clinical  response,  and  ER/PR  status 
of  the  13  patients  who  completed  the  protocol.  The  median  age  of  the  patient  group  was 
46  years  old  (range  31  to  70  years).  Of  these  13  patients,  eight  were  premenopausal  of 
whom  5  had  a  prior  history  of  taking  oral  contraceptives.  Five  patients  were 
postmenopausal  of  whom  2  were  on  hormone  replacement  therapy.  Histologic  diagnosis 
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was  obtained  by  ultrasound-guided  needle  core  biopsy  in  1 1  patients  and  by 
mammographic  stereotactic-guided  biopsy  in  2  patients.  Six  patients  were  diagnosed  with 
invasive  ductal  carcinoma  (IDC)  of  whom  4  also  had  ductal  carcinoma  in  situ  (DCIS). 
Two  patients  were  diagnosed  with  invasive  lobular  carcinoma  (ILC)  of  whom  1  also  had 
DCIS  and  the  other  also  had  both  DCIS  and  lobular  carcinoma  in  situ  (LCIS).  Four 
patients  were  diagnosed  with  IBC  of  whom  1  also  had  DCIS  and  another  also  had  LCIS. 
One  patient  was  diagnosed  with  both  IDC  and  ILC.  Eight  patients  were  ER  positive  and 
5  were  PR  positive.  Three  patients  had  additional  breast  lesions  with  LDo  <  2  cm. 

Because  these  additional  lesions  did  not  fit  the  voxel  selection  criteria  for  MRS,  they 
were  not  evaluated. 

The  median  LDo  was  3.8  cm  (range  2.1  to  9.5  cm).  The  median  LDf  was  2.7  cm 
(range  0  to  9.5  cm).  The  median  [tCho]o  was  4.5  mmol/kg  (range  0.9  to  8.5  mmol/kg)  for 
objective  responders  and  1.4  mmol/kg  (range  0.8  to  2.3  mmol/kg)  for  nonresponders. 
Overall,  [tChojo  was  significantly  higher  for  objective  responders  than  for  nonresponders 
(Wilcoxon  rank-sum  test,  p  =  0.03).  The  median  [tCho]24  was  3.4  mmol/kg  (range  0  to 
6.1  mmol/kg)  for  objective  responders  and  1.7  mmol/kg  (range  0.8  to  3.4  mmol/kg)  for 
nonresponders.  The  median  [tCho]f  was  0.5  mmol/kg  (range  0  to  2.4  mmol/kg)  for 
objective  responders  and  1.4  mmol/kg  (range  0  to  3.0  mmol/kg)  for  nonresponders. 

Objective  Response 

Eight  of  13  patients  experienced  an  objective  response  with  diminished  lesion 
size.  The  median  %ALDf  was  -56%  (range  -35  to  -100%),  and  the  median  %A[tCho]24 


11 


was  -30%  (range  -16  to  -100%).  All  8  patients  who  were  objective  responders  had 
[tChojo  >  [tCho]24>  [tCho]f  (Table  1).  Figure  1  shows  [tCho]24  and  LDf  as  a  percent  of 
the  baseline  for  the  objective  responders. 

Figure  2  shows  MR  data  from  an  objective  responder.  At  baseline,  [tCho]  was  4.6 
mmol/kg  and  LD  was  4.0  cm.  Twenty-four  hours  after  the  first  dose  of  AC  was 
administered,  [tCho]  dropped  to  3.7  mmol/kg  whereas  LD  remained  at  4.0  cm.  At  that 
point,  %A[tCho]24  was  -20%  and  it  was  predicted  that  the  patient  would  have  an  objective 
response.  After  the  fourth  dose  of  AC  (day  64),  [tCho]  decreased  to  0.9  mmol/kg  and  LD 
was  1.7  cm.  Based  on  %ALDf,  which  was  -58%,  the  patient  was  classified  as  an  objective 
responder.  This  patient  had  palpable  lymphadenopathy.  Since  treatment  with  paclitaxel  is 
regarded  as  standard  therapy  for  node-positive  patients  (32),  this  patient  was  continued 
on  systemic  therapy  with  paclitaxel  prior  to  surgery.  Interestingly,  24  hours  after  the 
second  dose  of  paclitaxel,  [tCho]  increased  to  4. 1  mmol/kg  and  LD  remained  at  1 .7  cm. 
Although  not  shown  in  the  figure,  LD  after  the  fourth  dose  of  paclitaxel  remained  at  1 .7 
cm,  whereas  [tCho]  increased  to  4.8  mmol/kg. 

No  Response 

Five  of  13  patients  were  nonresponders.  The  median  %ALDf  was  -13%  (range  0 
to  -23%)  and  the  median  %A[tCho]24  was  13%  (range  0  to  50  %).  All  5  patients  who 
were  nonresponders  had  a  [tChojo  <  [tCho]24.  With  respect  to  [tChojo,  [tChojf  was  lower 
in  2  patients,  higher  in  2  patients,  and  the  same  in  1  patient  (Table  1).  Figure  3  shows 
[tCho]24  and  LDf  as  a  percent  of  the  baseline  for  nonresponders. 
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Figure  4  shows  MR  data  from  a  nonresponder.  At  baseline,  [tCho]  was  1 .4 
mmol/kg  and  LD  was  2.9  cm.  Twenty-four  hours  after  the  first  dose  of  AC  was 
administered,  [tCho]  increased  to  2.1  mmol/kg  whereas  LD  remained  at  2.9  cm.  At  that 
point,  %A[tCho]24  was  50%  and  it  was  predicted  that  the  patient  would  be  a 
nonresponder.  After  the  fourth  dose  of  AC  (day  65),  [tCho]  was  0.9  mmol/kg  and  LD  was 
2.7  cm.  Based  on  %ALDf  of  only  -7%,  this  patient  was  classified  as  a  nonresponder. 

Correlation  between  %A[tCho]24  and  %ALDf 

A  significant  positive  correlation  between  %A[tCho]24  and  %ALDf  was  found 
(R  =  0.79,  p  =  0.001).  Figure  5  shows  a  graph  of  the  correlation  between  %A[tCho]24  and 
%ALDf.  Different  symbols  were  used  to  indicate  the  group  into  which  each  patient  was 
categorized:  objective  responder  or  nonresponder,  based  on  the  RECIST  criteria.  All  5 
patients  who  were  nonresponders  (black  squares)  had  a  %A[tCho]24  >  0.  All  8  patients 
who  had  an  objective  response  (white  squares)  had  a  %A[tCho]24  <  0.  %A[tCho]24  was 
significantly  different  between  objective  responders  and  nonresponders  (Wilcoxon  rank- 
sum  test,  p  =  0.007). 

Discussion 

This  pilot  study  shows  that  *H  MRS  appears  to  predict  clinical  response  in 
patients  with  locally  advanced  breast  cancer  within  24  hours  of  receiving  the  first  dose  of 
primary  systemic  therapy.  These  results  suggest  that  the  addition  of  high  field  !H  MRS 
can  offer  a  substantial  advantage  over  using  MRI  alone  in  predicting  response  to  PST. 
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Overall,  patients  who  were  objective  responders  had  a  higher  [tChojo  when 
compared  to  patients  who  were  nonresponders.  All  patients  who  were  objective 
responders  had  a  lower  [tCho]24  as  compared  to  [tChojo,  and  in  lesions  with  measurable 
[tCho]24 ,  a  further  decrease  in  [tCho]  was  observed  after  the  fourth  dose  of  AC.  Perhaps 
the  immediate  decrease  in  [tCho]  within  24  hours  after  the  first  dose  of  AC  reflects 
inhibition  of  cellular  proliferation  and  the  acute  cytotoxic  effect  of  chemotherapy.  All 
patients  who  were  nonresponders  had  either  no  change  or  a  higher  [tCho]24  as  compared 
to  [tChojo.  This  may  be  related  to  the  fact  that  certain  breast  neoplasms  exhibit  de  novo 
resistance  to  chemotherapy  and  as  a  result,  PST  with  the  AC  regimen  would  have  no 
effect  on  the  proliferating  neoplastic  cells. 

It  is  expected  that  quantitative  ’H  MRS  will  not  only  improve  the  accuracy  of 
MRI  in  detecting  breast  cancer,  but  as  shown  here,  it  may  be  used  with  MRI  to  assess 
response  to  PST  early  in  the  course  of  treatment.  Similar  to  our  study,  one  group  of 
researchers  found  that  within  a  majority  of  patients  who  underwent  PST  for  breast  cancer, 
there  was  a  disappearance  or  a  decrease  in  tCho  signal  (21).  However,  this  earlier  study 
did  not  categorize  which  patients  had  a  change  in  tCho  signal  nor  did  it  describe  when  the 
change  in  signal  occurred.  In  a  study  using  *H  MRS  to  evaluate  treatment  response  in 
pediatric  gliomas,  the  ratio  of  tCho  signals  from  tumor  and  normal  white  matter  was 
measured  (23).  In  brain  tumors  that  responded  clinically,  there  was  a  decrease  in  the  tCho 
ratio,  whereas  in  patients  who  did  not  respond,  there  was  no  change  or  an  increase  in  the 
tCho  ratio. 
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The  precise  mechanism(s)  as  to  why  neoplastic  tissues  exhibit  elevated  tCho 
levels  still  remains  unanswered.  It  has  been  proposed  that  increased  phosphocholine,  the 
primary  metabolite  responsible  for  the  tCho  peak  in  neoplastic  tissue,  is  a  result  of 
increased  synthesis  of  membranes  by  replicating  cells.  Elevated  tCho  may  reflect  a 
change  in  the  balance  between  biosynthetic  and  catabolic  pathways  in  which  choline 
compounds  serve  as  both  precursors  and  catabolites  (33).  Other  researchers  have 
observed  that  the  predominant  intracellular  mechanism  responsible  for  the  augmented 
levels  of  phosphocholine  is  due  to  up-regulation  of  choline  kinase  (34). 

Although  no  PST  regimen  other  than  AC  was  considered  in  this  study,  it  is 
interesting  to  note  the  pattern  of  change  in  [tCho]  for  the  patient  who  continued  treatment 
with  4  doses  of  paclitaxel  prior  to  surgery.  Despite  the  fact  that  this  individual  was  an 
objective  responder  to  AC,  the  patient  did  not  respond  to  paclitaxel.  Between  the  last 
dose  of  AC  and  the  second  dose  of  paclitaxel,  the  patient’s  [tCho]  increased  by  355%, 
from  0.9  to  4. 1  mmol/kg,  whereas  the  lesion  size  did  not  change.  The  patient’s  lesion 
size,  even  after  the  fourth  dose  of  paclitaxel,  remained  at  1 .7  cm,  while  the  [tCho] 
increased  further  by  17%,  from  4.1  to  4.8  mmol/kg.  It  is  reasonable  to  assume  that  early 
changes  seen  in  lH  MRS  may  even  be  used  in  predicting  clinical  response  between 
different  regimens.  Although  the  findings  of  this  study  showed  that  %A[tCho]24  was  an 
early  marker  for  response  to  the  specific  cytotoxic  agent  AC,  further  work  is  needed  to 
determine  whether  early  prediction  is  possible  with  other  agents. 
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Among  all  objective  responders  in  our  study,  the  smallest  change  in  [tCho]  within 
24  hours  was  a  decrease  of  16%.  The  ability  to  predict  response  depends  on  having 
sufficient  accuracy  and  reproducibility  in  the  MRS  measurement,  and  thus  future  studies 
using  multiple  baseline  scans  are  needed  to  establish  the  variance  of  [tCho] 
measurements.  We  believe  the  accuracy  of  a  [tCho]  measurement  is  limited  primarily  by 
signal-to-noise  ratio  (SNR)  which  has  been  shown  to  increase  at  least  linearly  with 
magnetic  field  strength  (35, 36).  Thus,  the  4T  scanner  used  here  offered  a  significant 
advantage  over  1 .5T.  Further  studies  are  needed  to  determine  whether  1 ,5T  scanners  can 
provide  the  SNR  necessary  to  measure  [tCho]  accurately.  However,  newer  3T  clinical 
scanners  are  expected  to  perform  similarly  to  the  4T  research  scanner  used  here.  An 
additional  limitation  of  MRS  is  the  inability  to  predict  clinical  response  in  breast  cancers 
that  are  small,  diffuse,  or  irregularly  shaped.  Due  to  SNR  limitations,  [tCho]  in  voxels  <  1 
cm3  may  be  difficult  to  measure  accurately.  Fortunately,  the  minimum  voxel  size  used 
for  monitoring  response  to  PST  will  typically  be  greater  than  1  cm3,  since  tumors  must  be 
at  least  2  cm  in  LD  (Stage  II  or  higher)  to  be  eligible  for  PST  by  current  criteria.  Breast 
tumors  with  thin  or  linear  morphology  can  pose  a  problem  for  the  spectroscopist  planning 
the  voxel,  particularly  when  intense  lipid  signals  are  contained  in  the  voxel  (24).  For 
example,  data  from  one  patient  were  excluded  from  the  analysis  due  to  an  inability  to 
plan  a  suitable  voxel  for  [tCho]  measurements  at  baseline  and  following  PST.  At 
baseline,  the  MRI  revealed  a  diffuse,  linearly  shaped  mass.  Due  to  the  thin  conformation 
of  this  lesion,  it  was  not  possible  to  place  a  rectangular  voxel  containing  less  than  33%  fat 
by  volume,  as  required  by  the  quantification  method  (25). 
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The  observation  of  consistent  trends  in  the  %A[tCho]24  between  objective 
responders  and  nonresponders  suggests  that  in  vivo  'H  MRS,  with  MRI,  may  be  a 
sensitive  indicator  in  predicting  clinical  response  as  early  as  24  hours  after  the  first  dose 
of  PST.  Monitoring  patients  on  PST  with  MRS  appears  to  offer  a  means  to  detect  the 
presence  of  viable  and/or  proliferating  neoplastic  cells.  With  the  possibility  of  early 
prediction  comes  the  benefit  of  immediate  assessment  in  tailoring  an  effective  regimen 
for  individual  patients.  Perhaps  early  changes  seen  in  !H  MRS  scans  of  patients  who 
undergo  PST  may  even  serve  as  an  indicator  of  risk  for  recurrent  or  metastatic  breast 
cancer.  Furthermore,  the  opportunity  for  such  early  prediction  will  be  useful  to 
researchers  seeking  to  evaluate  new  PST  drugs  and  may  help  to  elucidate  the  mechanisms 
of  multi-drug  resistance  in  breast  cancer. 

These  promising  findings  were  obtained  from  a  small  group  of  patients  and  thus  a 
study  of  a  larger  prospective  patient  series  is  needed.  A  larger  study  may  also  provide 
information  about  the  possibility  of  using  MRS  to  predict  when  a  patient  is  destined  to 
have  a  pathologic  complete  response.  To  the  best  of  our  knowledge,  this  is  the  first  report 
in  which  quantitative  tCho  measurements  were  made  with  in  vivo  'H  MRS  for  the 
purpose  of  predicting  clinical  response  to  PST  in  breast  cancer. 
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Table  1 


Patient 

No 

LDo 

(cm) 

LDf 

(cm) 

[tCho]o 

(mmol/kg) 

[tCho]  24 
(mmol/kg) 

[tCho]  f 
(mmol/kg) 

Diagnosis 

Clinical 

Response 

ER/PR 

1 

4.0 

1.7 

4.6 

3.7 

1.0 

IDC/DCIS 

OR 

+/- 

2 

6.8 

3.3 

1.7 

0.9 

0.6 

ILC/DCIS 

OR 

+/- 

3 

9.3 

4.2 

7.8 

5.4 

2.4 

IBC 

OR 

-/+ 

4 

2.8 

0 

2.7 

0 

0 

IDC/DCIS 

OR 

+/+ 

5 

3.1 

0.8 

0.9 

0.7 

0 

IDC 

OR 

-/- 

6 

2.1 

1.0 

4.5 

3.0 

0.4 

ILC/LCIS/DCIS 

OR 

+/+ 

7  + 

2.3 

0.7 

4.5 

3.8 

0 

IDC/DCIS 

OR 

-/+ 

8 

4.8 

3.1 

8.5 

6.1 

1.2 

IDC/DCIS 

OR 

+/- 

9 

9.5 

9.5 

1.1 

1.1 

0.3 

1BC/DCIS 

NR 

-/- 

10 

2.9 

2.7 

1.4 

2.1 

1.4 

IDC 

NR 

-/- 

11  + 

oo 

cn 

3.3 

0.8 

0.8 

0 

IDC/ILC 

NR 

+/+ 

12 

2.5 

2.0 

1.5 

1.7 

2.5 

IBC/LCIS 

NR 

+/- 

13  + 

4.7 

3.6 

2.3 

3.4 

3.0 

IBC 

NR 

+/- 

LDo  =  lesion  size  at  baseline,  LDf  =  lesion  size  after  the  fourth  dose  of  AC,  [tCho]o  =  total  choline 
concentration  at  baseline,  [tCho]24  =  total  choline  concentration  within  24  hours  after  the  first  dose  of  AC, 
[tCho]f  =  total  choline  concentration  after  the  fourth  dose  of  AC,  IDC  =  invasive  ductal  carcinoma,  ILC  = 
invasive  lobular  carcinoma,  IBC  =  inflammatory  breast  cancer,  DCIS  =  ductal  carcinoma  in  situ,  LCIS  = 
lobular  carcinoma  in  situ,  OR  =  objective  response,  NR  =  no  response,  ER  =  estrogen  receptor,  PR  = 

progesterone  receptor,  +  =  positive,  -  =  negative,  t =  patients  with  additional  satellite  lesions  less  than  2  cm 
in  longest  diameter  that  did  not  fit  the  inclusion  criteria  for  MRS  measurements. 
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Figure  Captions 

Figure  1  -  The  effects  of  treatment  on  the  concentration  of  choline-containing 
compounds  ([tCho])  and  lesion  size  (LD)  in  patients  with  an  objective  response.  The 
dashed  line  represents  the  normalized  value  of  [tCho]  and  LD  at  baseline.  The  diamonds 
represent  [tCho]  at  24  hours  after  the  first  dose  of  AC  normalized  to  baseline  and 
expressed  as  a  percentage  (%[tCho]24)  and  the  circles  represent  LD  after  the  fourth  dose 
of  AC  normalized  to  baseline  and  expressed  as  a  percentage  (%LDf).  %[tCho]24  and 
%LDf  pairs  from  individual  subjects  are  indicated  by  connecting  lines.  Notice  that  all  the 
objective  responders  had  a  [tCho]24  that  was  less  than  [tCho]  at  baseline. 

Figure  2  -  MRI  and  MRS  of  a  patient  with  an  objective  response  to  AC.  The  patient  was 
a  43  year-old  female  who  was  diagnosed  with  invasive  ductal  carcinoma  with  positive 
lymph  nodes.  Images  of  the  right  breast  (sagittal  view)  were  acquired  at  7  minutes  post 
Gd-DTPA  injection  using  fat-suppressed  3D  FLASH  (13.5/4.1).  The  box  surrounding  the 
enhancing  lesion  depicts  the  MRS  voxel.  Corresponding  spectra  are  shown  to  the  right. 
The  spectral  peaks  labeled  1-3  in  (a)  arise  from  lipid,  tCho,  and  water,  respectively.  The 
lines  above  and  below  the  tCho  peak  represent  the  fitted  tCho  peak  and  the  residual  of  the 
fit,  respectively,  (a)  Two  days  prior  to  starting  AC,  tCho  concentration  ([tChojo)  =  4.6 
mmol/kg  and  lesion  size  (LDo)  =  4.0  cm.  (b)  Within  24  hrs  after  the  first  dose  of  AC, 
tCho  concentration  ([tCho]24)  =  3.7  mmol/kg  and  lesion  size  (LD)  =  4.0  cm.  (c)  After  the 
fourth  dose  of  AC,  tCho  concentration  ([tChojf)  =  0.9  mmol/kg  and  lesion  size  (LDf)  = 

1 .7  cm.  (d)  After  the  second  dose  of  paclitaxel,  [tCho]  =  4. 1  mmol/kg  and  LD  =  1.7  cm. 
Using  equation  4,  the  percent  change  in  tCho  concentration  (%A[tCho]24)  was  -20%, 
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predicting  an  objective  response  to  AC.  Using  equation  2,  the  percent  change  in  lesion 
size  (%ALDf)  was  -58%,  compatible  with  an  objective  response  to  AC. 

Figure  3  -  The  effects  of  treatment  on  the  concentration  of  choline-containing 
compounds  ([tCho])  and  lesion  size  (LD)  in  nonresponders.  The  dashed  line  represents 
the  normalized  value  of  [tCho]  and  LD  at  baseline.  The  diamonds  represent  [tCho]  at  24 
hours  after  the  first  dose  of  AC  normalized  to  baseline  and  expressed  as  a  percentage 
(%[tCho]24)  and  the  circles  represent  LD  after  the  fourth  dose  of  AC  normalized  to 
baseline  and  expressed  as  a  percentage  (%LDf).  %[tCho]24  and  %LDf  pairs  from 
individual  subjects  are  indicated  by  connecting  lines.  Notice  that  all  the  nonresponding 
patients  had  a  [tCho]24  which  was  greater  than  or  equal  to  [tCho]  at  baseline. 

Figure  4  -  MR1  and  MRS  of  a  nonresponding  patient.  The  patient  was  a  42  year-old 
female  who  was  diagnosed  with  invasive  ductal  carcinoma.  Images  of  the  right  breast 
(sagittal  view)  were  acquired  at  7  minutes  post  Gd-DTPA  injection  using  fat-suppressed 
3D  FLASH  (13.5/4.1).  The  box  surrounding  the  enhancing  lesion  depicts  the  MRS  voxel. 
Corresponding  spectra  are  shown  to  the  right.  The  spectral  peaks  labeled  1-3  in  (a)  arise 
from  lipid,  tCho,  and  water,  respectively.  The  lines  above  and  below  the  tCho  peak 
represent  the  fitted  tCho  peak  and  the  residual  of  the  fit,  respectively,  (a)  One  day  prior  to 
starting  AC,  tCho  concentration  ([tChojo)  =1.4  mmol/kg  and  lesion  size  (LDo)  =  2.9  cm. 
(b)  Within  24  hrs  after  the  first  dose  of  AC,  tCho  concentration  ([tCho]24)  =  2.1  mmol/kg 
and  lesion  size  (LD)  =  2.9  cm.  (c)  After  the  fourth  dose  of  AC,  tCho  concentration 
([tChojf)  =  0.9  mmol/kg  and  lesion  size  after  the  fourth  dose  of  AC  (LDf)  =  2.7  cm. 
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Using  equation  4,  the  percent  change  in  tCho  concentration  (%A[tCho]24)  was  50%, 
predicting  no  response  to  AC.  Using  equation  2,  the  percent  change  in  lesion  size 
(%ALDf)  was  -7%,  compatible  with  a  nonresponder  to  AC. 

Figure  5  -  Correlation  between  percent  change  in  tCho  concentration  within  24  hours  of 
the  first  dose  of  AC  (%A[tCho]24)  and  percent  change  in  lesion  size  after  the  fourth  dose 
of  AC  (%ALDf).  Each  solid  white  square  represents  one  of  the  8  objective  responding 
patients  and  each  solid  black  square  represents  one  of  the  5  nonresponding  patients. 
There  is  a  statistically  significant  positive  correlation  between  %A[tCho]24  and  %ALDf, 
(R  =  0.79,  p  =  0.001). 
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